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SUMMARY 


This  report  presents  the  results  of  a  study  initiated 
by  the  Colorado  River  Water  Quality  Office  of  the 
USBR  (Bureau  of  Reclamation)  to  identify  possible 
systematic  changes  in  salt  loading  to  the  Colorado 
River  and  its  major  tributaries.  A  systematic  change 
may  involve  long-term  alterations  of  natural 
sources  of  salt  or  the  impact  of  human  activities  on 
salt  loading  or  salt  transport  in  the  river.  Systematic 
changes  have  been  suggested  as  a  possible  cause 
of  the  recent  decreases  in  the  salt  concentration  at 
Imperial  Dam  in  the  Lower  Basin.  Reducing  the 
salt  load  in  the  river  is  the  goal  of  the  Colorado 
River  Water  Quality  Improvement  Program.  If 
systematic  changes  in  the  salt  loading  character- 
istics of  the  basin  have  occurred,  they  must  be 
considered  in  formulating  plans  to  achieve  salinity 
reduction. 

A  systematic  change  might  manifest  itself  as  a 
long-term  shift  in  the  concentration,  load  (absolute 
mass),  or  mass  fraction  (relative  mass)  of  one  or 
more  of  the  major  solutes  in  the  system.  Probable 
shifts  can  be  identified  using  trend  analysis  or  by 
comparison  of  mean  values  before  and  after  major 
alterations  within  the  basin.  Both  approaches 
were  used  in  this  study.  Regression  techniques 
were  employed  to  test  for  significant  trends  and 
Student's  t  test  was  used  to  assess  the  mean  vaiue 
changes.  Six  major  ions  were  selected  for  analysis: 
calcium,  magnesium,  sodium,  chloride,  sulfate, 
and  bicarbonate.  Together  these  account  for  more 
than  95  percent  of  the  total  solute  load  at  the  major 
gage  sites  in  the  basin.  Fifteen  sites,  representing 
the  major  divisions  of  the  basin,  were  studied. 
Data  for  these  sites  were  obtained  from  the  USGS 
(U.S.  Geological  Survey)  computerized  water  data 
base,  WATSTORE. 

To  isolate  potential  sources  of  significant  trends  or 
mean  value  changes,  these  changes  must  be 
traceable  from  downstream  to  upstream  gage 
sites.  Therefore,  consistent  data  sets  were  required 
for  all  stations.  The  first  step  in  the  study  was  to 
create  a  data  file  of  representative  monthly  average 
streamflow  and  solute  concentration  values  for 
the  selected  gage  sites.  These  computer  data  files 
are  being  maintained  for  future  research  on  the 
Colorado  River  system. 

The  data  generated  was  first  used  to  test  various 
methods  relating  solute  concentration  to  stream- 
flow.  This  was  necessary  so  changes  in  concen- 
tration caused  by  natural  flow  variations  could  be 
accounted  for  in  subsequent  regression  analyses. 
A  theoretical  residence  time  model  was  compared 


with  the  widely  used  empirical  power  model.  The 
two  models  produced  equally  good  predictions  for 
periods  prior  to  flow  regulation.  These  models  are 
available  for  updating  input  data  used  in  the  USBR 
Colorado  River  Simulation  System  models. 

Although  15  gage  sites  were  selected,  trend 
analyses  were  completed  for  only  the  1 2  sites  that 
had  adequate  records  of  water  quality  prior  to 
major  flow  regulation.  Several  different  regression 
models  were  used,  and  their  results  were  generally 
corroborative.  Significant  negative  (decreasing) 
trends  were  found  for  magnesium,  sodium,  and 
sulfate  at  approximately  75  percent  of  the  sites.  At 
the  Lees  Ferry  and  Grand  Canyon  gage  sites, 
concentrations  of  these  three  ions  declined  in 
proportion  to  their  electrical  charge  and  relative 
mass,  resulting  in  a  reduction  in  TDS  (total  dis- 
solved solids)  of  approximately  2  mg/L  per  year 
for  the  low  flow  period,  December-March,  from 
1926  to  1961. 

After  construction  of  the  Colorado  River  Storage 
Project  reservoirs  in  the  early  1 960s,  the  variability 
of  both  solute  concentrations  and  mass  fractions 
was  reduced  at  downstream  gaging  stations.  This 
reduction  was  generally  proportional  to  the  degree 
of  regulation.  At  some  sites,  trends  toward  further 
reduction  in  variability  developed  after  initial  reser- 
voir filling.  Several  other  activities  of  man,  including 
mining,  urbanization,  and  salinity  control  practices, 
may  have  affected  water  quality  in  the  basin 
during  this  period.  Impacts  decreased  with  distance 
downstream,  but  were  often  identified  at  more 
than  one  gaging  station.  However,  the  major 
reservoirs,  Lakes  Powell  and  Mead,  concealed  all 
impacts  that  had  occurred  in  the  Upper  Basin. 

No  significant  trends  in  the  load  for  any  solute 
have  been  identified  for  the  station  below  Hoover 
Dam  since  1 965;  however,  the  concentration  of  all 
solutes  except  chloride  decreased  there  during  at 
least  one  season. 

Additional  research  is  recommended  to  determine 
the  possibility  and  probable  magnitude  of  con- 
tinuing impacts  of  major  reservoirs  on  down- 
stream salinity.  Also,  because  of  the  relative 
chemical  stability  at  the  gage  site  below  Hoover 
Dam,  the  changes  observed  downstream  at 
Imperial  Dam  may  have  resulted  from  changing 
conditions  in  the  Lower  Basin  area.  Therefore,  an 
extension  of  the  present  study,  with  additional 
Lower  Basin  gage  sites  included,  is  recommended. 
This  extension  should  also  update  the  monthly 
data  base  and  data  base  generation  procedures, 
and  limit  statistical  analyses  to  those  recommended 
in  the  present  work. 


INTRODUCTION 

Description  of  the  Basin 

The  Colorado  River  Basin  includes  632  000  km2 
in  the  Western  United  States  and  northern  Mexico 
(fig.  1).  The  average  unregulated  flow  of  the 
Colorado  River  below  Lees  Ferry,  Arizona,  is 
16.0-18.5  krnVa,  which  is  small  compared  with 
that  of  other  North  American  rivers  with  similar 
size  basins.  However,  the  river  is  an  important 
source  of  water  for  more  than  12  million  people 
and  approximately  1  million  ha  of  irrigated  agri- 
cultural land  [1].* 

The  headwaters  of  the  Colorado  and  its  major 
tributaries,  the  Green  and  San  Juan  Rivers,  lie  in 
the  high  peaks  of  the  Rocky  Mountains,  where  the 
annual  precipitation  is  normally  between  100  and 
1 50  cm.  Most  of  its  course,  however,  crosses  the 
semiarid  Colorado  Plateau  and  the  Sonoran  Desert, 
where  the  average  annual  precipitation  is  only 
6  cm  [2].  Many  of  the  geologic  formations  in  this 
part  of  the  basin  are  of  marine  origin  and  contain 
sodium  chloride  (halite)  and  calcium  sulfate 
(gypsum)  salts.  Natural  springs  and  man-made 
wells  intercept  saline  ground  waters  associated 
with  these  formations  and  discharge  into  the  river 
system.  Soils  in  much  of  the  basin  have  developed 
residually  on  gypsum-bearing  shales.  Irrigation 
water  applied  to  this  land  promotes  weathering 
and  dissolution  of  salts  from  the  soil  and  under- 
lying shales.  As  a  result  this  water  returns  to  the 
river  with  a  greater  salt  load  than  was  diverted  [3]. 
Irrigation  also  increases  dissolved  mineral  con- 
centrations in  the  river  by  depleting  the  streamf  low 
volume.  As  a  result,  TDS  (total  dissolved  solids) 
increases  from  approximately  50  mg/L  at  the 
headwaters  to  800  mg/L(1 977-1 981  average)[2] 
at  Imperial  Dam,  the  final  diversion  point  on  the 
Colorado  River  in  the  United  States. 

Water  resource  development  in  the  basin  began 
when  the  pioneers  settled  there  in  the  1 860's.  By 
the  1920's,  when  the  USGS  (U.S.  Geological 
Survey)  began  monitoring  surface  water  quality, 
much  of  the  present  irrigation  development  was 
already  in  place,  and  several  reservoirs  and  trans- 
basin  diversions  had  been  constructed  in  the 
upper  reaches  of  the  basin.  The  first  major  multi- 
purpose reservoir,  Lake  Mead,  began  filling  in 
1 935.  Development  in  the  basin  upstream  of  Lake 
Mead  continued  gradually,  though  mainstem  flow 
was  essentially  unregulated  until  the  early  1 960's. 
Between  1962  and  1966,  the  storage  capacity 
above  Lees  Ferry  increased  from  3  to  45  km3, 


*  Numbers  in  brackets  refer  to  entries  in  the  bibliography. 


mostly  as  a  result  of  the  construction  of  the 
Colorado  River  Storage  Project  reservoirs:  Navajo 
Reservoir  on  the  San  Juan  River  (began  filling  in 
November  1962),  Blue  Mesa  Reservoir  on  the 
Gunnison  River  (November  1 965),  and  Lake  Powell 
on  the  Colorado  River(March  1963).  Reservoir  and 
transbasin  diversion  development  in  the  basin 
upstream  of  Lake  Mead  is  outlined  in  table  1  of 
appendixes  B  through  P. 


Colorado  River  Water  Quality  Improvement 
Program 

During  the  period  of  accelerated  development,  the 
USGS  initiated  an  appraisal  of  the  water  supply  in 
the  basin  [3].  Its  purpose  was  to  determine  whether 
development  would  be  limited  by  "legal,  physical, 
and  economic  factors.  "Much  of  the  study  focused 
on  the  impacts  of  human  activities  on  salinity.  It 
concluded  that,  by  1957,  approximately  half  the 
average  annual  TDS  concentration  at  Lees  Ferry 
was  caused  by  domestic,  industrial,  and  agri- 
cultural activities  within  the  basin  and  transbasin 
diversions. 

By  1 970,  salinity  in  the  lower  Colorado  River  was 
recognized  as  a  basinwide  problem  [4].  TDS  stan- 
dards for  lower  basin  gage  locations  were  adopted 
by  the  basin  States  in  response  to  amendments  to 
the  Federal  Water  Pollution  Control  Act  of  1972. 
These  standards  were  set  at  the  1972  average 
concentrations,  thus  establishing  a  policy  of  non- 
degradation.  To  offset  the  effects  of  future  water 
development,  the  Colorado  River  Basin  Salinity 
Control  Act  of  1974  authorized  the  USBR  to 
construct  4  salinity  control  projects  and  investigate 
the  feasibility  of  constructing  12  others.  The 
location  of  active  projects  upstream  of  Hoover 
Dam  is  shown  in  figure  1 . 

At  the  same  time,  computational  methods  were 
developed  to  predict  the  effects  of  planned  water 
resources  development.  Early  projections  indi- 
cated that  salinity  increases  observed  at  Imperial 
Dam  between  1949  and  1970  would  continue, 
with  TDS  concentration  eventually  reaching 
1200  mg/L  [4,  5,  6].  However,  contrary  to  these 
projections,  salinity  has  actually  decreased  since 
1970  (see  fig.  2). 

It  is  not  clear  whether  this  decrease  in  salinity  is 
part  of  a  long-term  cycle  or  is  indicative  of  a 
permanent  change  somewhere  in  the  system. 
Other  investigators  [7,  8]  have  argued  that  new 
conditions  created  by  impoundment  of  large 
reservoirs  have  permanently  altered  chemical 
processes  in  the  river  system  and  may  be  a  factor 


A  Solinity  Control   Project  Sites 
O  Goge  Stotion : 

I.     Colorado  R.  neor   Cameo ,  Co 

Gunnison  R.   near  Grand   Junction,  Co. 
Dolores  R    near  Cisco,  Ut. 
Colorado  R.  near  Cisco,  Ut. 
Green  R.  neor  Green  River,  Wy. 
Green  R    neor  Greendole,  Ut. 
Yampa  R.  neor  Maybell,  Co. 
Duschene  R    near   Randlett,  Ut. 
White  R    neor   Watson,  Ut. 
Green  R.  at  Green   River,  Ut. 
San  Rafael   R.  near  Green    River,  Ut. 
neor    Bluff,  Ut. 
at   Lee's   Ferry,  Az. 
below   Grond  Conyon.Az. 
below  Hoover  Dam.Az. 


2. 

3 

4 

5. 

6. 

7. 

8. 

9. 

10. 

II. 

12. 

13. 

14. 

15. 


San  Juon  R 
Colorado  R. 
Colorado  R. 
Colorado  R 


Las  Vegas 


MEXICO 

Figure  1 .  —  Map  of  salinity  control  project  sites  and  gage  stations  in  the  Colorado  River  Basin. 
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in  decreasing  salinity.  Other  possible  factors  are 
manmade  or  natural  alterations  in  the  salt  load 
producing  areas  of  the  basin.  At  a  given  point  in 
the  system,  a  systematic  change  might  manifest 
itself  as  a  shift  in  the  solute  concentration- 
streamflow  relationship  or  in  the  absolute  or 
relative  solute  load. 


Research  Objectives 

As  a  first  step  in  evaluating  the  possibility  of  a 
systematic  change,  the  present  study  was  initiated 
with  the  following  objectives: 

1 .  Create  a  data  base  that  includes  average  month- 
ly streamflows  and  major  solute  concentra- 
tions at  selected  gaging  stations  in  the 
Colorado  River  Basin 

2.  Develop  a  theoretical  model  for  the  relation- 
ship between  streamflow  and  individual 
solute  concentrations 

3.  Identify  significant  long-term  trends  in  solute 
concentration,  load,  and  mass  fraction  for 


the  period  before  the  construction  of  major 
upstream  dams 

4.  Identify  significant  changes  in  the  solute 
concentration,  load,  and  mass  fraction  after 
the  accelerated  reservoir  development  of  the 
early  1960's 

Fifteen  stations  were  selected  for  analysis  (table  1 ). 
These  are  all  principal  stations  used  to  define 
concentration-streamflow  relationships  in  the 
USBR's  Colorado  River  Simulation  System.  Except 
for  the  station  on  the  Duchesne  River  and  those 
below  Lake  Mead  and  Flaming  Gorge  Reservoir,  all 
had  adequate  water  quality  records  for  the  planned 
statistical  analyses. 


This  report  describes  the  procedures  used  to 
develop  the  data  base  and  presents  the  results  of 
the  statistical  analyses.  It  includes  an  evaluation  of 
the  selected  statistical  procedures  based  on  the 
ease  of  interpreting  their  results.  A  cursory  inter- 
pretation of  the  results  and  their  implications  for 
the  salinity  control  program  is  also  given,  though  a 
more  complete  analysis  is  recommended.  No 
attempt  is  made  to  reconcile  the  recently  observed 
TDS  at  Imperial  Dam  with  earlier  projections. 


Table  1 .  —  Colorado  River  Basin  gaging  stations  selected  for  analysis 


USGS 

number 


Location 


Period  of 
record1 


90955 
91525 
91800 
91805 
92170 
92345 
92510 
93020 
93065 
93150 


Colorado  River  near  Cameo,  Colorado 
Gunnison  River  near  Grand  Junction,  Colorado 
Dolores  River  near  Cisco,  Utah 
Colorado  River  near  Cisco,  Utah 


Green  River  near  Green  River,  Wyoming 
Green  River  near  Greendale,  Utah 
Yampa  River  near  Maybell,  Colorado 
Duchesne  River  near  Randlett,  Utah 
White  River  near  Watson,  Utah 
Green  River  at  Green  River,  Utah 


93285         San  Rafael  River  near  Green  River, 
93795         San  Juan  River  near  Bluff.  Utah 


Utah 


12/33- 

10/31- 

5/51- 

10/28- 


6/79 
8/82 
9/82 
4/83 


10/51- 
10/56- 
1 2/50- 

1/51- 
1 2/50- 

8/28- 


5/82 
6/83 
8/82 
3/83 
7/79 
4/83 


2/47- 
fi/?fl- 


5/83 
10/81 


.93800        Colorado  River  at  Lees  Ferry,  Arizona 

94025 

94215 


Colorado  River  near  Grand  Canyon,  Arizona 
Colorado  River  below  Hoover  Dam,  Arizona 


5/27- 


11/25- 
10/40- 


6/82 
"3/74 
4/83 


Includes  all  data  necessary  for  generation  of  monthly  concentrations  as  described  in  this  report. 


DATA  BASE  DEVELOPMENT 

A  series  of  computer  procedures  was  developed  to 
generate  representative  monthly  values  from 
existing  data.  These  procedures  were  also  used  to 
check  for  errors  in  the  data  base  and  to  estimate 
missing  values.  A  complete  guide  to  the  use  of  the 
component  programs  is  given  in  appendix  A. 


attempted;  but  the  sum  of  six  ions,  necessary  to 
compute  mass  fraction,  was  estimated  from 
regression  on  evaporation  residue  TDS.  The 
periodic  chemical  analyses,  including  missing  value 
estimates,  were  then  combined  with  daily  values 
of  streamflow  and  conductance  to  estimate  solute 
concentrations  for  days  samples  were  not  taken. 


WATSTORE  Retrieval 

All  data  used  in  this  project  were  retrieved  directly 
from  the  USGS  WATSTORE  data  base  by  an 
automatic  communication  link  between  the  USBR 
CYBER  and  USGS  IBM  computers.  These  data 
i nclude  mean  daily  streamflows,  daily  specifi c 
conductances,  and  periodic  chemical  analyses  of 
the  following:  calcium,  magnesium,  sodium,  chlo- 
ride, sulfate,  bicarbonate,  alkalinity,  TDS  (residue 
on  evaporation  at  180  °C),  and  specific  conduct- 
ance. Chemical  analyses  are  reported  for  time- 
weighted  composite  samples  until  September  1 969 
and  for  monthly  discrete  samples  thereafter. 


Generation  of  Monthly  Mean  Values 

Composite  sample  analyses  were  assigned  to 
each  day  of  the  composite  period;  discrete  sample 
analyses  were  assigned  to  the  day  of  collection. 
Substitutions  for  missing  daily  values  weregener- 
ated  by  interpolation  in  a  hierarchical  process  in 
which  streamflow  was  estimated  first,  followed  by 
conductance,  and  finally  solute  concentrations. 
This  sequence  was  employed  because  streamflow 
was  used  to  weight  the  interpolation  of  conduc- 
tance, and  conductance,  in  turn,  was  used  to 
weight  the  interpolation  of  concentration.  The 
general  formula  used  was: 


Initial  Analysis 


**=(k 


X>o+ 


k-k0 

Kn  —  K, 


iwknxkn  - 


wkxko)\  wk-  (1) 


Potential  errors  in  the  periodic  chemical  data  were 
identified  by  testing  the  differences  between  the 
observed  and  predicted  values  of:  (1)  total  solute 
concentration,  for  which  the  predicted  value  was 
the  sum  of  constituents  and  the  observed  value 
was  the  residue  on  evaporation  (increased  to 
account  for  lost  bicarbonate);  (2)  specific  conduc- 
tance, for  which  the  predicted  value  was  computed 
from  equivalent  ionic  conductances;  and  (3)  net 
charge,  for  which  the  predicted  value  was  assumed 
to  be  zero.  Obvious  errors  were  checked  and 
corrected,  as  described  in  appendix  A. 

Possible  errors  in  daily  values  were  identified  from 
outliers  in  the  log-log  regression  of  conductance 
on  streamflow.  Typographical  errors  in  the 
WATSTORE  data  were  confirmed  and  corrected  by 
comparing  them  with  values  in  USGS  Water 
Supply  Papers.  Other  obviously  erroneous  values 
were  set  to  zero  and  treated  as  missing  values  in 
subsequent  procedures.  Changes  made  to 
WATSTORE  data  are  documented  in  table  2  of 
appendixes  B  through  P. 

Single  missing  solute  concentrations  in  individual 
chemical  analyses  were  estimated  assuming  a  net 
charge  of  zero.  When  two  or  more  values  were 
missing  from  an  analysis,   no  estimates  were 


where:     X      =  Streamflow,  conductance,  or 

solute  concentration 
W     =  A  weighting  factor 
k       -  Julian  date  within  a  period  of 

missing  values 
k0     -  Last  Julian  date  with  observed 

data  before  the  period  of  missing 

values 
kn     =  First  Julian  date  with  observed 

data  after  the  period  of  missing 

values 


Weighting  factors  in  equation 
follows: 


1   are  defined  as 


1 .  Streamflow:  Wk  -  1 ,  for  all  k 

2.  Conductance:  Wk  =  1  /Qk  ,  where  B  is  the 
slope  of  the  log-log  regression  of  daily  con- 
ductance on  daily  streamflow  (Q)  for  the 
entire  period  of  unregulated  flow.  Wk  -  1 ,  if  B 
is  not  significantly  different  from  0  (because 
of  major  flow  regulation). 

3.  Solute  concentrations:  Wk  =  1  /Lk,  where/,  is 
specific  conductance  on  day  k. 


For  streamflow  interpolation  equation  1  is  linear 
with  respect  to  time. 

Conductance  and  concentration  are  transformed 
by  the  weighting  factors  into  parameters  (ak  and  bk 
in  equations  2  and  3)  that  are  expected  to  be 
relatively  constant  over  short  periods.  For  inter- 
polation of  conductance,  equation  1  becomes  a 
linear  interpolation,  with  respect  to  time,  of  the 
parameter  ak  in  the  equation: 


Lk  -  ak  Q* 


(2) 


Lane  [9]  found  this  power  relationship  to  be 
generally  applicable  to  natural  streams  in  the 
Western  United  States.  Conductance  interpolation 
was  not  attempted  when  the  missing  value  period 
included  days  of  no  flow. 

The  interpolation  of  solute  concentration  is  a 
linear  interpolation  of  the  parameter  bk  in  the 
equation: 


Ck  -  bt  L, 


(3) 


Over  small  ranges,  this  linear  relationship  provides 
a  good  estimate  of  concentration.  After  inter- 
polation, solute  concentrations  were  adjusted  to 
account  for  differences  between  the  instantaneous 
specific  conductance  (Ls)  at  the  time  of  sample 
collection  and  the  mean  daily  conductance  (LD).  if 
the  latter  value  was  available. 


Ck  -  Ck 


(4) 


The  effect  of  the  length  of  the  missing  value  period 
on  the  accuracy  of  solute  concentration  estimates 
was  tested  using  data  from  the  discrete  sampling 
period  for  the  Colorado  R  iver  at  Lees  Ferry,  Arizona. 
Increasing  the  maximum  period  from  45  to  60 
days  resulted  in  change  of  only  3  percent  in  the 
average  monthly  concentrations  for  the  sum  of  6 
ions.  In  addition,  the  number  of  months  with 
inadequate  data  for  computing  monthly  means,  as 
defined  below,  was  reduced  by  21  percent.  To 
maintain  representative  data  on  the  maximum 
number  of  months,  the  interpolation  interval  was 
set  at  60  days,  which  allows  interpolation  between 
any  two  samples  collected  in  consecutive  calendar 
months.  Testing  the  accuracy  of  streamflow  inter- 
polation was  considered  unnecessary  because  of 
the  infrequent  occurrence  of  gaps  in  these  data. 

Following  interpolation,  monthly  mean  values  of 
streamflow  and  flow-weighted  concentrations 
were  computed  from  daily  values.  When  measured 
or  estimated  values  were  unavailable  for  more 
than  25  percent  of  the  days  in  a  month,  analysis  of 
that  month  was  discontinued.  Therefore,  the 
resultant  monthly  values  are  intended  to  preserve 
the  maximum  amount  of  information  from  existing 
data  without  masking  long-term  trends  or  including 
unreasonable  estimates.  Monthly  values  are  listed 
in  table  3  of  appendixes  B  through  P. 


Analysis  of  Seasonality 

Prior  to  statistical  analysis,  the  generated  monthly 
data  were  tested  for  seasonal  variability.  Normal- 
ized deviations  (A)  from  the  annual  mean  were 
computed  for  each  month  using  the  formula: 


This  adjustment  was  made  only  for  dates  with 
continuously  recorded  specific  conductance.  It 
was  made  to  correct  potentially  nonrepresentative 
instantaneous  sample  concentrations  to  average 
conditions  on  the  sample  date. 

Interpolations  were  limited  to  intervals  of  20  days 
for  streamflow  and  conductance,  and  60  days  for 
solute  concentrations.  The  effect  of  the  length  of 
the  missing  value  period  on  the  accuracy  of 
specific  conductance  estimates  was  tested  using 
data  from  the  Colorado  River  at  Cisco,  Utah. 
Artificial  gaps  of  varying  length  were  created  in  the 
conductance  record,  and  generated  estimates  were 
compared  with  measured  values.  Standard  errors 
were  6.8  and  10.2  percent  of  the  observed  mean 
for  7-day  and  20-day  intervals,  respectively. 
Although  increasing  the  interpolation  interval 
decreased  the  estimation  accuracy,  a  10  percent 
standard  error  was  not  considered  unreasonable. 


A  = 


Xy,m        Xy 


(5) 


where:    X    =  streamflow,  solute  concentration, 

or  load 
X    =  the  annual  mean  of  each 

parameter 
y     =  year 
m  =  month 


The  average  monthly  values  for  A  were  then 
plotted  to  determine  the  seasonal  pattern.  Figure  3 
shows  the  seasonal  variation  in  streamflow  for  the 
Colorado  River  at  Cisco,  Utah.  A  peak  runoff  period 
in  May  and  June  is  clearly  indicated.  April  and  July 
are  transition  months,  and  the  remaining  months 
appear  to  constitute  a  baseflow  period.  However, 
chemical  variation,  shown  in  figure  4,  indicates  a 
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Figure  3.  —  Average  monthly  streamflow  variation,  Colorado  River  at  Cisco,  Utah,  1928-1982. 
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Figure  4.  —  Average  monthly  variation  in  TDS  load,  Colorado  River  at  Cisco,  Utah,  1928-1982. 


difference  between  the  August-November  and 
December-March  periods.  Higher  concentrations 
in  the  former  are  mainly  caused  by  calcium  and 
sulfate,  which  are  predominant  in  irrigation  return 
flow.  Other  gage  sites  gave  similar  results,  unless 


streamflow  was  regulated  by  an  upstream  reser- 
voir. Therefore,  the  annual  cycle  was  divided  into 
three  seasons:  ( 1 )  basef  low,  December-March,  (2) 
peak  runoff,  May-June,  and  (3)  return  flow, 
August-November. 


CONCENTRATION-STREAMFLOW 
RELATIONSHIPS 

In  unregulated  streams  salinity  generally  varies 
inversely  with  streamf  low.  Possible  causes  of  this 
inverse  relationship  and  a  review  of  several  mathe- 
matical models  proposed  to  describe  this  relation- 
ship are  listed  by  Lane  [9].  In  evaluating  the 
beneficial  effects  of  salinity  control  projects, 
concentration-streamflow  models  are  useful  to 
correct  the  salinity  for  the  effect  of  the  natural 
variability  of  streamf  low  before  and  after  construc- 
tion of  the  project.  In  this  report,  these  models  are 
used  in  multiple  regression  analyses  to  remove 
predictable,  short-term  fluctuations  in  concen- 
tration, thus  allowing  more  precise  estimates  of 
long-term  changes  in  concentration  and  load.  This 
chapter  develops  a  theoretical  model  with  two 
parameters  that  have  physical  meaning  and 
compares  the  results  of  this  model  with  those 
obtained  from  the  power  equation. 

The  power  equation  is  an  empirical  model  com- 
monly used  to  describe  the  concentration-stream- 
flow  relationship  [9].  This  equation  (with  solute 
subscripts  omitted)  is: 


4.  The  solubility  of  any  solute  is  constant;  i.e., 
the  solubility  is  independent  of  the  chemical 
composition  of  both  the  dissolving  solids  and 
the  aqueous  solution. 

The  differential  equation  (with  solute  subscripts 
omitted)  for  dissolution  is: 


QdC  =  -k(C'-CsJdA' 


(7) 


where:     Q    =  Streamf  low  (m3/s) 

C"  =  Solute  concentration  after  Q  has 
contacted  A',  the  area  of  solid  sur- 
face (g/m3) 
Cs  =  Solute  solubility  (g/m3) 
k     =  Mass  transfer  coefficient  (m/s) 
A'  =  Solid-liquid  contact  area  (m2) 

By  assumption  1,  the  boundary  condition  is: 

CfA' =  0)  =  0  (8) 

By  assumptions  3  and  4,  the  solution  to  equation  7 
is: 


C  =  Cs[1-exp(-r)] 


(9) 


C  =  bjQ 


b2 


(6) 


where  t,  the  dimensionless  residence  time,  is 
defined  by: 


where:     C    =  Concentration 
Q    -  Streamflow 
bj,  b2    =  Empirical  parameters  estimated 
by  regression 

The  parameter  b2  ranges  from  zero,  where  con- 
centration is  independent  of  streamflow,  to  minus 
one.  where  load  is  independent  of  streamflow. 
Equation  6  yields  unrealistically  high  estimates  of 
concentrations  at  low  flows,  because  C  approaches 
infinity  as  Q  approaches  zero. 

Because  the  power  model  fails  at  low  flows  and 
because  the  parameter  b1  and  b2  have  no  physical 
meaning,  a  theoretical  residence  time  model  is 
derived  based  on  the  assumptions  listed  below. 
Assumptions  2  and  4  may  represent  serious 
oversimplifications  of  the  physical  system. 

1 .  The  solute  concentration  in  rain  or  snowmelt 
is  negligible. 

2.  The  contact  area  of  dissolvable  solids  is 
independent  of  the  streamflow. 

3.  The  mass  transfer  coefficient  for  the  dis- 
solution of  solids  is  constant. 


_  kA 
Q 


(10) 


and  A  is  the  total  solid-liquid  contact  area  upstream 
of  the  gaging  station.  Equation  9  is  of  the  form: 


C  =  £,  [1-exp  H>2/Q)] 


(11! 


where  b1  and  b2  are  parameters  to  be  estimated  by 
regression,  based  on  observed  values  of  concen- 
tration and  streamflow.  The  parameter  b1  repre- 
sents the  solute  solubility,  and  b2  represents  the 
product  of  the  mass  transfer  coefficient  and  the 
total  solid-liquid  contact  area.  The  residence  time 
equation  yields  realistic  values  of  concentrations 
at  zero  streamf  low  (infinite  residence  time),  where 
the  solute  concentration  approaches  solubility, 
and  at  infinite  streamflow  (zero  residence  time), 
where  the  concentration  approaches  zero  (see 
fig.  5). 

To  test  the  effectiveness  of  the  power  and  resi- 
dence time  equations  in  describing  the  relationship 
between  monthly  values  of  concentration  and 
streamflow,  regressions  were  conducted  for  6 
solutes  and  their  sum.  Regressions  were  run  for 
each  of  the  three  seasons  at  12  different  gaging 
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Figure  5.  —  Residence  time  equation  on  log-log  scale.  (The 
power  equation  would  yield  a  straight  line  with  the  slope 
equal  to  the  power  term  b2  in  equation  6.) 


stations  on  the  Colorado  River  (see  tables  4  and  5 
in  appendixes  B  through  0).  Nonlinear  regressions 
of  equations  6  and  1 1  used  the  SPSS  (Statistical 
Package  for  the  Social  Sciences)  subprogram 
NONLINEAR  [11,1 2]  with  the  Marquardt  algorithm 
[13].  Values  only  from  periods  of  minimum  flow 
regulation  upstream  of  the  gaging  stations  are 
included  in  the  regression  analyses.  By  dividing 
the  year  into  three  seasons,  assumptions  2  and  4 
for  the  residence  time  model  implythatA,  andb2ir\ 
equation  1 1  are  constant  only  for  one  season. 
Although  the  solute  solubility  parameter  (b,  in 
equation  1 1 )  may  be  independent  of  the  season, 
seasonal  values  were  estimated  to  maintain  the 
same  degrees  of  freedom  used  in  the  power 
regression. 

Because  precision  in  the  chemical  analyses  is 
generally  proportional  to  the  magnitude  of  the 
solute  concentrations,  a  weighting  factor  of  Cm"2  is 
used  in  the  nonlinear  regressions.  With  this 
weighting  factor,  the  RMS  (root-mean-square) 
error,  or  standard  deviation  of  the  residuals,  is 
dimensionless  and  becomes  the  fractional  standard 
deviation.  For  the  weighted  regression,  the  RMS 
error  is  defined  as: 


"-,       /  Cj.obs  ~  Cj, 


RMS  error  = 


pre. 


Cj. 


iobs 


n-k 


0.5 


(12) 


where:        Cjobs  -  Observed  monthly  value  of 
concentration 


Cj.pred  ~  Monthly  value  of 

concentration  estimated  from 
streamflow  by  the  power  or 
residence  time  equation 
n  =  Number  of  observations 
k  =  Number  of  estimated 
parameters  (k  -  2  for 
equations  6  and  1 1) 
n-k  =  degrees  of  freedom 

An  RMS  error  of  0. 1 0  indicates  (assuming  random 
normally  distributed  residuals  and  greater  than  30 
degrees  of  freedom)  an  approximately  68-percent 
probability  that  the  estimated  and  observed  values 
of  concentrations  differ  by  less  than  1 0  percent. 

Although  complete  model  evaluation  requires 
examination  of  the  residuals  for  any  systematic 
bias  in  estimating  concentration,  residual  exami- 
nation is  omitted  from  this  report  because  of  the 
large  number  of  regressions.  A  preliminary  com- 
parison of  the  accuracy  of  the  models  can  be  made 
from  the  RMS  errors  summarized  in  table  2.  The 
missing  values  in  table  2  result  from  difficulties  in 
computing  that  regression.  Based  on  the  RMS 
errors,  the  power  equation  is  more  accurate  than 
the  residence  time  equation  for  the  peak  flow 
months  of  May  and  June.  For  the  other  two 
seasons  the  RMS  errors  differ  little.  Neither  equa- 
tion is  accurate  for  the  Dolores  River,  near  Cisco, 
Utah,  where  the  RMS  errors  are  generally  greater 
than  0.20.  Both  equations  are  used  in  the  multiple 
regression  analyses  in  the  next  section  to  estimate 
long-term  changes  in  concentration. 

The  regression  results  of  the  power  equation  (see 
table  4  in  appendixes  B  through  0)  indicate  that 
for  several  solutes  at  several  locations,  the  stream- 
flow  has  only  a  minor  effect  on  the  concentration. 
Based  on  the  estimated  values  of  the  power  term 
(b2  in  equation  6)  that  are  greater  than  -0.20,  the 
bicarbonate  concentration  is  relatively  indepen- 
dent of  the  streamflow  for  all  stations  (see  table  3). 
By  the  same  criteria,  the  calcium  concentration 
varies  little  with  the  streamflow  at  7  of  the  12 
stations. 


LONG-TERM  CHANGES  FOR 

PERIODS  OF  MINOR  FLOW 

REGULATION 

The  object  of  this  section  is  to  determine  whether 
the  monthly  values  of  concentration,  load,  or  mass 
fraction  of  any  of  the  six  major  ions  (calcium, 
magnesium,  and  sodium  cations;  and  chloride, 
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Table  2.  —  Root-mean-square  errors  (weighted  by  Cm~2)  of  two  equations  describing  the  concentration-streamflow 

relationship 

RMS  error  of  Cm  =  brQjb2  /  RMS  error  of  Cm  =  6,[1  -  exp(-b2/Qm)] 

Sum  of 
Period  Season  Calcium         Magnesium         Sodium  Chloride  Sulfate         Bicarbonate         six  ions 


Colorado  River  near  Cameo,  Colorado 

1933-45  Dec-Mar         0.057/0.052  0  068/0  062  0  064/0  058  0.066/0.064  0.066/0  062  0.054/0  049  0.053/0.049 

May-Jun.         0  127/0.130  0.241/  0.078/  0.085/0  087  0.070/0  070  0.098/0.098  0.052/0.053 

Aug. -Nov.        0.045/0.066  0  239/0  226  0.188/  0.063/0  068  0.067/0.070  0.039/0.045  0.042/0.058 

Gunnison  River  near  Grand  Junction,  Colorado 

1931-64  Dec-Mar.        0.107/0107  0111/0.110  0  142/0.142   0.165/0164  0.121/0.120  0.082/0.082  0103/0.103 

May-Jun.         0.108/0185  0151/0206  0180/0233  0  303/0  327  0.125/0165  0.090/0.116  0.094/0.152 
Aug-Nov.        0.121/0.125  0172/0156  0.176/0179  0187/0205  0  145/0.141    0.097/0.084  0.128/0.129 

Dolores  River  near  Cisco,  Utah 

1951-82  Dec -Mar.        0.324/0  324  0  348/0  343   0  379/0  371    0  393/0  382  0.334/0.329  0.823/  0.337/0.332 

May-Jun.         0.198/0253  0  225/0  306  0.061/  0.244/  0.200/0.280  0.281/0.289  0181/0.268 

Aug-Nov         0.237/0  270  0  184/0  199  0.344/0  342  0.396/0  392  0  288/0  304  0.947/  0.221/0.236 

Colorado  River  near  Cisco,  Utah 

1928-46  Dec -Mar.        0.084/0.087  0.103/0.105  0.109/0110  0.137/0.137  0.115/0.118  0.065/0.067  0.085/0.087 

May-Jun.         0.092/0.117  0.146/0.163  0.107/0.105  0.122/0.119  0.127/0.145  0.126/0.121    0.087/0.115 
Aug-Nov.        0.097/0.117  0  090/0.105   0  114/0.117  0.245/0.254  0.114/0.127  0.080/0  066  0.077/0.096 


Green  Basin 


Green  River  near  Green  River,  Wyoming 


1951-62  Dec-Mar.         0.095/0.095  0.155/0.152   0.149/0.154  0.140/0.149  0.198/0.199  0.089/0.092  0.137/0.139 

May-Jun.         0.135/0.162  0.217/0.240  0  138/0.164  0.197/0.223  0.170/0.192  0.139/0.151    0.139/0.172 

Aug-Nov.        0.182/0.182  0.158/0.157  0.150/0.150  0.179/0.179  0.194/0.192  0.106/0.106  0.139/0.138 

Yampa  River  near  Maybell,  Colorado 

1950-82  Dec-Mar.  /0.112  0201/0.196  0.183/0.165  0.256/0.246  0.291/  0.102/0.101    0.140/0.124 

May-Jun.         0.241/  0.392/0.396  0.336/0.341    0.376/0.388  0.376/  0.265/  0.283/ 

Aug. -Nov.        0.124/0.122  0.194/0.182   0.215/0.200  0.301/0.295  0.223/0.219  0.129/0.126  0.145/0.145 

White  River  near  Watson,  Utah 

1950-79  Dec-Mar.        0.108/0.108  0.114/  0.249/  0.400/  0.131/  0.091/  0.155/ 

May-Jun.         0.122/0.131  0.148/0.171    0.252/0  261    0.340/0.328  0.173/0.196  0.128/0.139  0.129/0.150 

Aug.-Nov.        0.103/  0.131/0.134  0.229/  0.384/  0.140/  0.091/  0  117/ 

Green  River  at  Green  River,  Utah 

1928-61  Dec-Mar.        0.079/0.086  0.122/0.117  0.104/0.105  0.110/0.130  0.1 13/0.1 1 1   0.073/0.077  0.079/0.082 

May-Jun.         0.112/  0.205/  0.206/0.214  0.151/0.164  0.214/0.220  0.132/  0.144/0.145 

Aug.-Nov.        0.110/  0.184/0.180  0.166/0.168  0.174/0.183  0.180/0.180  0.193/0.092  0.118/0.120 
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Table  2.  —Root-mean-square  errors  (weighted  by  Cm2)  of  two  equations  describing  the  concentration-streamflow 

relationship  —  Continued 


Sum  of 
Period  Season  Calcium         Magnesium         Sodium  Chloride  Sulfate         Bicarbonate         six  ions 

San  Rafael  River  near  Green  River,  Utah 

1947-65  Dec-Mar         0  105/  0  136/  0  150/  0.136/0.161    0.134/0.148  0.146/  0  115/ 

May-Jun.         0.134/0.237  0.159/0291    0.216/0289  0.183/0286  0.191/0.280  0.060/0062  0.152/0.263 
Aug. -Nov         0.175/0.202  0  320/0.308  0.423/0.436  0.251/0.291    0  208/0.218  0.204/0.200  0.194/0.203 


San  Juan  Basin 


San  Juan  River  near  Bluff,  Utah 


1929-40          Dec -Mar         0073/0.101    0  108/0  128  0  153/0  165  0.137/0  156  0.116/0.138  0  070/0  075  0088/0.110 

May-Jun          0.100/0.120  0.139/0.141    0.208/0.247  0.324/0  143  0.150/0  182  0  098/0  103  0.107/0  138 

Aug. -Nov         0  149/0  185  0  186/0  213  0.259/0.301    0  167/0  213  0  211/0  255  0.125/0.131  0.177/0.218 

San  Juan  River  near  Bluff,  Utah 

1929-61           Dec. -Mar.        0.078/0.104  0  129/0.139  0.155/0  175  0.172/0.183  0.110/0.132  0.073/0076  0.090/0.114 

May-Jun          0.136/0.157  0.219/0.232  0.210/0.236  0.192/0212  0.175/0.202  0.142/  0.146/0.172 

Aug.-Nov.        0143/0.166  0.232/0.235  0  222/0.264  0.188/0223  0189/0228  0.131/0.130  0.157/0.196 

Lower  Colorado  Basin 

Colorado  River  at  Lees  Ferry,  Arizona 

1927-41           Dec-Mar.        0.055/0  075  0  036/0.052  0.056/0  074  0.056/0  066  0.052/0.069  0.040/0.055  0.040/0  061 

May-Jun.         0.028/0.026  0  082/0.094  0.208/0  201    0.051/0.049  0.074/0.064  0.059/0.062  0.059/0.058 

Aug.-Nov         0.191/             0.138/0.151    0.170/0194  0.145/0.157  0211/0222  0.066/0.070  0.165/0.181 

Colorado  River  at  Lees  Ferry,  Arizona 

1927-61           Dec-Mar.        0.071/0079  0124/0.123  0.063/0.068  0059/0061  0  079/0  084  0.071/0.073  0.051/0.058 

May-Jun          0.129/0.143   0.164/0.179  0.159/0.173  0.144/0  156  0.128/0  147  0  133/0.140  0.105/0.127 

Aug -Nov.        0.152/0.164  0.167/0.176  0  152/0.155  0  150/0  149  0  167/0.179  0.096/0.095  0  126/0.142 

Colorado  River  near  Grand  Canyon,  Arizona 

1925-41           Dec-Mar.        0.063/0.071    0.086/0.086  0.07C/0.078  0  064/0  064  0.075/0  080  0  059/0.064  0.057/0.064 

May-Jun          0.111/0128  0127/0.140  0.133/0149  0.104/0  108  0  147/0  166  0107/0.112  0.106/0.127 

Aug.-Nov         0149/             0.122/0.155  0.129/0164  0.088/0.107  0169/0200  0072/0.090  0.126/0.163 

Colorado  River  near  Grand  Canyon,  Arizona 

1925-61           Dec-Mar.        0.057/0  064  0  100/0  099  0  078/0  086  0  144/0.147  0.080/0.084  0.064/0.066  0.056/0  062 

May-Jun          0.143/              0.153/0.166  0.139/0.111    0.098/0.104  0.144/0.164  0.167/0.168  0.125/0.140 

Aug -Nov.         0.142/0.153  0  138/0.154  0.129/0.147  0.093/0.105  0.164/0.182  0.098/0.100  0.118/0.139 
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Table  3.  —  Identification  of  stations,  seasons,  and  solutes  where  solute  concentration  (C)  or  load(L) 

is  relatively  independent  of  streamflow 


Period 


Season         Calcium     Magnesium      Sodium         Chloride 


Sum  of 
Sulfate      Bicarbonate      six  ions 


Upper  Colorado  Basin 


Colorado  River  near  Cameo,  Colorado 

1933-45  Dec-Mar. 
May-Jun. 
Aug. -Nov. 

Gunnison  River  near  Grand  Junction,  Colorado 


1931-64 


Dec. -Mar. 
May-Jun. 
Aug. -Nov. 


Dolores  River  near  Cisco,  Utah 

1951-82       Dec-Mar.  C 

May-Jun. 
Aug. -Nov. 

Colorado  River  near  Cisco,  Utah 

1928-46  Dec-Mar. 
May-Jun. 
Aug. -Nov. 


Green  River  near  Green  River,  Wyoming 

1951-62       Dec-Mar.  C  C 

May-Jun.  C 

Aug. -Nov. 

Yampa  River  near  Maybell,  Colorado 


1950-82       Dec-Mar.  C 

May-Jun.  C 

Aug. -Nov.  C 


White  River  near  Watson,  Utah 

C 
C 

Green  River  at  Green  River,  Utah 

1928-61 


1950-79  Dec-Mar. 
May-Jun. 
Aug. -Nov. 


Dec-Mar. 
May-Jun. 
Aug. -Nov 


Green  Basin 


c 
c 


c 
c 
c 


c 
c 


c 
c 
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Table  3.  —  Identification  of  stations,  seasons,  and  solutes  where  solute  concentration  (C)  or  load(L)  is 

relatively  independent  of  streamflow  —  Continued 

Sum  of 
Period  Season         Calcium     Magnesium      Sodium         Chloride         Sulfate      Bicarbonate      six  ions 


San  Rafael  River  near  Green  River,  Utah 


1 947-65 

Dec. -Mar. 
May-Jun. 

Aug. -Nov. 

C 

C 

San  Juan 

River  near  Bluff, 

Utah 

1 929-40 

Dec. -Mar. 

C 

May-Jun. 

C 

C 

Aug. -Nov. 

C 

C 

San  Juan 

River  near  Bluff, 

Utah 

1929-61 

Dec. -Mar. 

May-Jun. 

C 

Aug. -Nov. 

C 

c 
c 

San  Juan  Basin 


Lower  Colorado  Basin 


Colorado  River  at  Lees  Ferry,  Arizona 

1927-41        Dec-Mar.  C 

May-Jun. 
Aug. -Nov.  C 

Colorado  River  at  Lees  Ferry,  Arizona 

1927-61  Dec-Mar. 
May-Jun. 
Aug. -Nov.  C 

Colorado  River  near  Grand  Canyon,  Arizona 

1925-41  Dec -Mar. 
May-Jun. 
Aug. -Nov.  C 

Colorado  River  near  Grand  Canyon,  Arizona 
1925-61        Dec-Mar.  C 

C 


C 
C 


Dec-Mar. 
May-Jun. 
Aug. -Nov. 


C 
C 
C 


C 
C 
C 


C 
C 


C 
C 
C 


C 
C 


In  the  power  equation: 

Cm  =  b1Qmb2 

the  solute  concentration  is  identified  as  being  relatively  independent  of  streamflow  when  the  estimated 
value  of  b2  is  greater  than  -0.20  and  the  weighted  RMS  error  is  less  than  0.50.  For  solute  load, 
independence  is  identified  for  values  of  b2  less  than  -0.80. 
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sulfate,  and  bicarbonate  anions)  had  changed 
significantly  before  construction  of  the  major  up- 
stream dams.  Monthly  values  for  the  solute  con- 
centration are  the  flow-weighted  means  of  daily 
concentrations  estimated  by  the  procedure  de- 
scribed under  DATA  BASE  DEVELOPMENT.  The 
monthly  solute  load  is  computed  by: 


J,m  =  CKmQm/1000 


(13) 


where: 


Ji.n 
C,r 

Qrr 


=  Monthly  value  of  load  of  solute  / 

in  kg/s 
=  Monthly  value  of  concentration 

of  solute  /in  mg/L 
=  Monthly  value  of  streamflow  in 

m3/s  computed  as  the  mean  of 

daily  measurements  (see  DATA 

BASE  DEVELOPMENT) 
=  Solute  (Ca*2,  Mg*2,  Na2,  CI, 

S04  2,  HCO3,  or  tne  sum  of  these 

six  ions) 


The  monthly  solute  mass  fraction  Wim  expressed 
as  a  percent  of  each  of  the  six  major  ions  is 
computed  by: 


r 

W,m  =  -r *  100  percent 

^  7.m 


(14) 


where:     C7m    =  sum  of  the  concentrations  of 
the  six  major  ions  (refer  to 
DATA  BASE  DEVELOPMENT) 

Because  concentration  and  load  are  generally 
strong  functions,  and  mass  fraction  may  be  a  weak 
function  of  streamflow,  multiple  regression  analy- 
ses are  used  to  determine  the  existence  of  long- 
term  trends  in  these  parameters.  In  determining 
the  existence  of  long-term  changes  in  salinity 
during  periods  of  minor  flow  regulation,  the  fol- 
lowing equations  (with  solute  subscripts  omitted) 
can  be  used  in  the  regression  analyses: 


Jm  =  b7+  b2Qm  +  b3t 

Cm  =  b7+  b2Qm  +  b3t 

Wm  =  b,+  b2Qm  +  b3t 

Cm-biQmb>tb* 

Wm  =  blQmb2tb3 


(15) 
(16) 
(17) 
(18) 
(19) 


Cm  =  b,Qmb2exp(b3t) 


Wm  =  b1Qmb2exp(b3t) 


6;[1-exp(-tf>2  +  /M/QJ] 

(weighting  factor  =  Cm 


(20) 


(21) 


(22) 


where:     bh  b2,  b3 


-  Parameters  estimated  by 

regression.  Different 

values  of  these 

parameters  are  estimated 

for  each  equation. 
=  Year  in  the  twentieth 

century  (e.g.,  for  the  year 

1 950,  t  =  50). 


For  evaluating  the  statistical  significance  of  long- 
term  changes  in  the  water  quality  parameters,  the 
regression  paramenter  b3  in  each  of  the  above 
eight  equations  is  tested  against  the  null  hypo- 
thesis that  b3  is  zero. 

Equations  15,  16,  and  17  are  standard,  un- 
weighted, empirical  linear  equations  for  use  in 
multiple  linear  regression.  The  computed  statis- 
tical significance  of  b2  in  equation  15  is  ignored 
because  of  the  intrinsic  correlation  between 
streamflow  and  solute  load,  which  is  computed  as 
the  product  of  monthly  values  of  streamflow  and 
concentration.  The  estimated  long-term  linear 
trend  for  the  dependent  variable  y  (Jm,  Cm,  or  Wm) 
in  equations  15,  16,  and  17  is: 


5y 
8t 


=  b. 


(23) 


where  the  subscript  Qm  indicates  that  the  derivative 
is  evaluated  at  a  constant  value  of  Qm. 

Equations  18  and  19  are  unweighted  empirical 
power  equations,  whose  logarithms  are  linear  and 
solvable  by  multiple  linear  regression.  Power 
equations  are  often  used  for  parameters  with  log- 
normal  distributions.  The  power  equation  for  con- 
centration (equation  18)  is  used  because  of  the 
accuracy  of  the  power  equation  (equation  6)  for 
describing  the  concentration-streamf  low  relation- 
ship. The  estimated  rate  of  change  of  the  de- 
pendent variable  y  (Cm  or  Wm)  in  equations  1 8  and 
19  is: 


= b3y/t 


(24) 


15 


and  the  average  long-term  trend  can  be  estimated 
by  evaluating  equation  24  at  the  mean  values  of  y 
and  t.  The  average  long-term  trend  for  load  asso- 
ciated with  the  trend  in  concentration  can  be 
estimated  bv: 


^    )       =b3QmCm/(W00t)        (25) 


8t 


where  b%  is  estimated  by  regression  with  equation 
1 8  and  Qm,  Cm,  and  t  are  the  mean  values  of  Qm, 
Cm,  and  t. 

Equation  20  was  also  derived  from  the  empirical 
power  relationship  by  allowing  br  inequation  5  to 
vary  with  time.  With  a  linear  change  in  bu  equation 
6  may  be  written  as: 


Cm  =  b,(1  +b3t)Qmb2 


(26) 


which  can  be  approximated  by  equation  20  for 
absolute  values  of  b3t  less  than  1 .0.  Multiple  linear 
regression  can  be  used  with  the  logarithms  of 
equations  20  and  21  to  estimate  the  parameters 
bj,b2,  and  b3.  The  rate  of  change  of  the  dependent 
variable  y  (Cm  or  Wm)  in  equation  20  and  21  is: 


=  b3y 


(27) 


and  the  average  long-term  trend  can  be  estimated 
by  evaluating  equation  29  at  the  mean  values  of 
Cm  and  Qm.  The  average  long-term  trend  for  the 
load  is: 


&4 
8t 


=  b3(b1-CJ/W00 


(30) 


Equations  29  and  30  are  valid  approximations  for 
the  average  long-term  trends  only  wheji  the 
estimated  value  of  b1  is  greater  than  Cm.  In 
equations  25, 28,  and  30,  the  estimated  long-term 
trends  for  the  load  associated  with  trends  in 
concentration  are  intended  to  be  independent  of 
the  long-term  changes  in  streamflow.  Changes  in 
the  load  associated  with  changes  in  streamflow 
are  estimated  by: 


&4 
8t 


1000     8t 


(31 


In  regression  equation  1 5  for  load  the  parameter 
b3  measures  both  effects,  and  thus  includes  the 
additional  changes  in  load  due  to  transmountain 
diversions  or  other  diversions  with  no  return  flow. 
The  long-term  trend  for  the  load  associated  with  a 
linear  change  in  only  concentration  can  be  esti- 
mated from  the  regression  results  of  equation  16 
by: 


and  the  average  long-term  trend  can  be  estimated 
by  evaluating  equation  27  at  the  mean  value  of  y. 
The  average  long-term  trend  for  load  associated 
with  the  trend  in  concentration  can  be  estimated 
by: 


8Jn 

8t 


^P  =b?QmCm/1000 


(28) 


where  b3  is  the  regression  parameter  in  equation 
20. 

Equation  22  is  the  residence  time  equation  with  an 
additional  regression  parameter,  b3,  for  long-term 
changes  in  residence  time  caused  by  changes  in 
the  solid-liquid  mass  transfer  coefficient  or  in  the 
contact  area.  Equation  22  is  not  linear  with  respect 
tothe  parameters/);,  i^'31"1*^;}  and  must  be  solved 
by  weighted,  nonlinear  regression.  The  rate  of 
change  of  Cm  is: 


8Cn 
8t 


=  b3(b1-CJ/Qr 


(29) 


8t 


=  b3Qm/'\000 


(32) 


Multiple  linear  regressions  of  equations  15 
through  21  used  FORTRAN  subroutines  BECOVM 
and  RLMUL  of  the  IMSL  (International  Mathe- 
matical and  Statistical  Libraries,  Inc.)  Library  [14]. 
Nonlinear  regression  of  equation  22  employed 
the  SPSS  subprogram  NONLINEAR  with  the 
Marquardt  algorithm.  Complete  results  of  the  re- 
gression analyses  of  equations  1 5  through  20  are 
contained  in  tables  6  through  9  of  appendixes  B 
through  0. 

Tables  4  through  7  summarize  the  long-term 
trends  for  the  concentration  and  associated  long- 
term  trends  for  the  load  estimated  by  regression  of 
concentration  equations  1 6, 1 8,  20,  and  22  where 
the  parameter  b3  is  significantly  different  from 
zero  at  the  95-percent  confidence  level.  Table  8 
lists  the  statistically  significant  trends  with  the 
highest  confidence  levels  from  tables  4  through  7. 
The  four  equations  yielded  similar  results  with 
several  cases  where  the  estimated  long-term 
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trends  are  identical  to  two  decimal  places.  Separate 
estimates  of  b1  for  each  season  would  result  in  a 
better  fit  to  the  data  and  in  more  precise  estimates 
of  long-term  trends  by  equations  30  and  31, 
where  the  difference  between  b,  and  Cm  is  evalu- 
ated. In  several  cases  in  table  7,  Cm  is  larger  than 
bj  (see  equations  30  and  31 ),  causing  a  reversal  in 
the  sign  of  the  estimated  long-term  trends.  In 
these  cases,  only  the  sign  of  b3  in  equation  22  is 
listed. 

Tables  9,  10,  and  1 1  contain  the  estimated  long- 
term  trends  for  mass  fraction  where  significant 
values  of  b3  are  indicated  in  equations  17,  19,  and 
21.  The  three  empirical  equations  estimated 
almost  identical  long-term  trends  in  mass  fractions. 
Table  1 2  contains  all  trends  that  are  indicated  to  be 
significant  by  any  of  the  three  mass  fraction 
regressions. 

Where  no  trend  is  indicated  for  an  entry  in  tables  4 
through  1 2,  an  actual  trend  may  exist,  but  the  data 
do  not  justify  rejection  of  the  null  hypothesis  of 
zero  long-term  trend.  Only  a  few  (4  to  1 5)  months 
of  data  are  available  for  the  Cameo  station,  which 
was  included  in  the  trend  analysis  because  of  its 
importance  as  a  major  gaging  station  on  the 
Colorado  River.  Generally,  for  each  row  in  tables  4 
through  8,  the  sum  of  the  significant  trends  in 
concentration  for  the  six  individual  ions  (columns 
3  through  8)  is  approximately  equal  to  the  signifi- 
cant trend  estimated  for  the  sum  of  the  six  ions 
(column  9).  For  each  row  of  mass  fractions  in 
tables  9  through  12,  the  sum  of  the  significant 
trends  of  the  individual  ions  is  approximately  zero. 

Long  term  declines  in  the  TDS  (sum  of  six  ions), 
corrected  for  the  variation  in  streamflow  by  multiple 
regression,  are  evident  at  8  of  the  12  gaging 
stations  (see  table  8).  The  Dolores  River  near  Cisco 
and  the  Yampa  River  near  Maybell  are  the  only 
stations  indicated  to  be  experiencing  an  increase 
in  salinity. 

At  eight  gaging  stations  the  regression  results 
indicate  a  negative  trend  in  sulfate  ion  concen- 
tration and  mass  fraction  (see  tables  8  and  12). 
Declines  in  the  concentrations  of  magnesium  and 
sodium  ions  are  evident  at  seven  stations  (not 
including  Cameo).  For  the  base  flow  season  of 
December  to  March,  in  the  Colorado  River  at  Lees 
Ferry  and  near  the  Grand  Canyon,  the  concen- 
trations of  the  sulfate,  magnesium,  and  sodium 
ions  have  declined  in  stoichiometric  proportions, 
causing  a  decline  in  TDS  of  approximately  2mg/L 
per  year.  Figures  6  through  9  and  1 0  through  1 3 
illustrate  the  declines  in  magnesium  and  sulfate 
concentrations  during  December,  January,  Feb- 


ruary, and  March  in  the  Colorado  River  near  the 
Grand  Canyon.  The  /-variable  is  transformed  to 
accommodate  a  two-dimensional  figure.  The  in- 
variable in  figures  6  through  8  and  1 0  through  1 2 
is  the  ion  concentration  with  a  linear  or  power 
correction  for  the  natural  variability  of  streamflow. 
The  regression  lines  are  the  results  of  the  empirical 
multiple  linear  regressions  (equations  1 5, 1 7,  and 
19).  The  /-variable  in  figures  9  and  13  is  the 
estimated  value  of  the  product  of  the  solid-liquid 
mass  transfer  coefficient  (k)  and  the  total  solid- 
liquid  contact  area  (A)  upstream  of  the  gaging 
station  (see  equations  8  and  9).  The  regression 
lines  are  the  result  of  nonlinear  regression  with 
the  residence  time  equation  (equation  22) 


IMPACT  OF  FLOW  REGULATION 

When  an  impoundment  is  created,  upstream 
waters  are  mixed  in  the  reservoir  and  quality 
fluctuations  diminish.  Water  released  downstream 
during  periods  of  natural  low  flow  decreases  the 
normally  high  salt  concentrations.  In  addition, 
chemical  changes  caused  by  the  dissolution  and 
precipitation  of  mineral  salts  within  the  reservoir 
may  change  the  character  of  the  release  water  [8]. 
At  locations  such  as  Lees  Ferry,  where  flow  is 
almost  completely  regulated,  concentration  is  no 
longer  related  to  flow.  At  other  sites,  where  flow 
has  become  partially  regulated,  the  concentration- 
flow  relationship  may  be  weakened  or  changed 
altogether.  The  impact  of  regulation  must  then  be 
evaluated  independent  of  streamflow.  Therefore, 
two  statistical  analyses  were  used  to  test  for 
changes  in  water  quality  as  a  function  of  time 
alone. 

In  both  analyses,  the  monthly  data  were  divided 
into  two  periods:  (1)  before  1963,  and  (2)  after 
1 965.  These  represent  the  times  before  and  after 
the  filling  of  the  Colorado  River  Storage  Project 
reservoirs.  Solute  loads,  mass  fractions,  and  flow- 
weighted  concentrations  were  computed  for  the 
three  seasons:  December-March,  May-June,  and 
August-November.  Student's  t  test  was  then  used 
to  determine  significant  differences  between  mean 
seasonal  values.  Trends  following  regulation  were 
tested  by  linear  regression,  with  time  as  the  sole 
independent  variable.  All  stations  were  treated 
alike,  even  if  the  flow  there  did  not  become 
regulated  during  the  1963-65  period.  The  results 
could  then  be  compared  between  stations  and 
regulation  impact  could  be  separated  from  other 
basinwide  influences.  Complete  results  of  these 
tests  are  displayed  in  tables  10  and  1 1  of 
appendixes  B  through  P. 
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Figure  6.  Linear  trend  of  magnesium  ion  concentration 
with  a  linear  correction  for  the  variability  of  stream- 
flow  for  the  Colorado  River  near  the  Grand  Canyon, 
December  through  March,  1925-61. 


Year    in   the  20th  Century 

Figure  8.  Semilog  trend  of  magnesium  ion  concentration 
with  a  power  correction  fur  the  variability  of  stream- 
flow  for  the  Colorado  River  near  the  Grand  Canyon, 
December  through  March,  1925-61 . 


Year    in   the  20th  Century 

Figure  7.  Log-log  trend  of  magnesium  ion  concentration 
with  a  power  correction  for  the  variability  of  stream- 
flow  for  the  Colorado  River  near  the  Grand  Canyon, 
December  through  March,  1925-61. 


Year    in  the  20th  Cr>ti.r» 

Figure  9.  Magnesium  ion  concentration  corrected  for 
streamflow  by  the  residence  time  equation.  The 
ordinate  is  the  product  of  the  mass  transfer  coef- 
ficient (k)  and  the  solid-liquid  contact  area  (-4)  up- 
stream of  the  gaging  station  on  the  Colorado  River 
near  the  Grand  Canyon,  December  through  March, 
1925-61. 
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Figure  1 0.  Linear  trend  of  sulfate  ion  concentration  with 
a  linear  correction  for  the  variability  of  streamf  low  for 
the  Colorado  River  near  the  Grand  Canyon,  December 
through  March,  1925-61. 


Year   in  the  20th  Century 

Figure  12.  Semilog  trend  of  sulfate  ion  concentration 
with  a  power  correction  for  the  variability  of  stream- 
flow  for  the  Colorado  River  near  the  Grand  Canyon, 
December  through  March,  1925-61. 


Year   in  the  20th  Century 

Figure  1 1 .  Log-log  trend  of  sulfate  ion  concentration 
with  a  power  correction  for  the  variability  of  stream- 
flow  for  the  Colorado  River  near  the  Grand  Canyon, 
December  through  March,  1925-61 . 


Year   in  the  20th  Century 

Figure  13.  Sulfate  ion  concentration  corrected  for 
streamflow  by  the  residence  time  equation.  The 
ordinate  is  the  product  of  the  mass  transfer  coef- 
ficient (k)  and  the  solid-liquid  contact  area  {A)  up- 
stream of  the  gaging  station  on  the  Colorado  River 
near  the  Grand  Canyon,  December  through  March, 
1925-61. 
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Changes  in  Mean  Values 

Significant  changes  in  concentration  and  mass 
fraction  for  selected  stations  are  shown  by  season 
in  tables  1 3  through  1 5.  From  1 963  through  1 965, 
major  reservoirs  were  created  on  the  Gunnison 
and  San  Juan  Rivers,  on  the  Green  River  above 
Greendale,  Wyoming,  and  on  the  Colorado  River 
above  Lees  Ferry,  Arizona.  The  obvious  impact  has 
been  a  general  decrease  in  solute  concentrations 
during  months  of  historically  low  flow  and  high 
concentration  (August-March)  and  an  increase 
during  what  had  been  peak  flow,  low  salinity 
months  (May-June).  This  results  from  mixing 
waters  of  different  quality  in  the  reservoir,  thereby 
reducing  fluctuations.  This  effect  was  also  con- 
firmed by  comparing  the  seasonal  variation  in 
pre-1963  with  post-1965  mean  values.  An 
example  showing  sulfate  concentration  and  mass 
fraction  is  given  in  table  16. 

For  all  stations  below  new  impoundments,  the 
range  of  seasonal  mean  values  is  greatly  reduced 
following  flow  regulation.  This  is  most  obvious  for 
stations  where  flow  is  completely  regulated.  At 
Lees  Ferry,  the  range  of  seasonal  average  sulfate 
concentrations  was  103-447  mg/L  before  regu- 
lation and  225-257  mg/L  after  regulation.  The 
range  of  sulfate  mass  fraction  was  reduced  from 
27-43  percent  to  38-39  percent.  A  similar  impact 
can  be  seen  for  the  Greendale  station.  Smaller,  but 
still  significant,  impacts  are  identified  for  the 
partially  regulated  stations:  Grand  Junction;  Cisco; 
Green  River,  Utah;  and  Bluff.  At  unregulated 
stations  (Maybell  and  Watson),  the  ranges  in- 
creased or  remained  approximately  constant. 

Associated  with  the  reduction  in  variability  is  a 
general  shift  in  relative  mass  from  sulfate  to 
bicarbonate  in  August  through  March,  and  vice 
versa  in  May  and  June.  The  reduction  in  sulfate 
mass  fraction  can  be  attributed  to  the  mixing  in  the 
reservoirs.  Sulfate  is  more  prevalent  in  base  and 
return  flows  because  its  major  source  is  dissolution 
of  gypsum  in  ground  water.  Therefore,  its  domi- 
nance during  these  periods  is  reduced  by  mixing 
with  the  lower  concentration  spring  runoff  water. 
There  was  also  a  relative  decrease  in  calcium 
during  May  and  June.  This  may  result  from  calcite 
precipitation  within  the  reservoirs  [8,  1 5]. 

Some  of  the  impacts  of  impoundment  can  be 
traced  downstream  to  the  stations  at  Cisco  and 
Green  River,  Utah.  Generally,  their  magnitude  and 
significance  are  somewhat  reduced  as  intervening, 
unregulated  inflow  exerts  a  moderating  effect. 
This  is  not  true,  however,  for  solute  loads  at  Cisco. 


Figure  14  displays  changes  in  streamflow  and 
anion  loads  as  they  occur  spatially  throughout  the 
basin.  Note  the  decreased  sulfate  load  at  Cisco 
during  the  base  and  return  flow  periods.  There 
does  not  seem  to  be  any  upstream  source  for 
either  reduction.  Streamflow  has  not  changed 
significantly  in  the  return  flow  season  and  has 
actually  increased  during  base  flow.  The  sulfate 
load  has,  therefore,  decreased  independently  of 
flow,  somewhere  on  the  Colorado  River  between 
Cameo  and  Cisco.  The  most  likely  cause  of  this  is  a 
reduction  in  loading  from  agricultural  sources  in 
the  Grand  Valley  near  Grand  Junction,  Colorado. 
Here,  also,  gypsum  dissolution  is  a  major  source  of 
low  flow  sulfate  load.  The  reduction  in  sulfate 
loading  could  be  the  result  of  two  activities:  (1 )  the 
salinity  control  practices  initiated  by  the  Soil 
Conservation  Service  in  the  mid-1 970's,  and  (2) 
the  urbanization  of  approximately  10  percent  of 
the  irrigated  area.  These  activities  involved  lining 
laterals  and  farm  ditches  and  reducing  irrigation 
applications.  Both  activities  could  reduce  subsur- 
face return  flow  volume  and  salt  pickup. 

The  changes  identified  for  the  station  at  Green 
River,  Utah,  cannot  be  attributed  entirely  to  up- 
stream impoundments.  Here,  also,  water  quality  is 
affected  by  intervening  inflows,  particularly  from 
the  Yampa  and  White  Rivers.  Although  neither  of 
these  tributaries  is  regulated,  both  exhibit  signifi- 
cant changes  in  water  quality.  The  concentrations 
and  relative  masses  of  sulfate  and  magnesium 
increased  during  all  three  seasons  at  the  Yampa 
station.  This  may  have  resulted  from  surface  coal 
mine  development  within  the  basin.  Other  investi- 
gators [6]  have  documented  sulfate  and  mag- 
nesium loading  from  coal  spoils  in  this  area. 

At  the  White  River  station,  chloride  and,  to  a  lesser 
extent,  sodium  decreased  in  both  concentration 
and  relative  mass;  these  decreases  were  most 
prominent  during  low  flow  seasons.  These  changes 
might  have  been  caused  by  the  plugging  of  Meeker 
Well,  upstream  of  Meeker,  Colorado,  in  1 968.  This 
well  was  a  point  source  of  sodium-chloride  con- 
tamination. It  flowed  at  a  regular  rate,  which  was 
greatly  reduced  for  1 1  years  after  it  was  plugged. 
This  flow  reduction  may  explain  the  greater  impact 
during  low  flow,  when  the  well  would  have 
contributed  a  greater  relative  water  volume  to  the 
river.  These  impacts  can  be  traced  downstream  to 
the  Green  River,  Utah  station. 

Downstream  of  the  stations  at  Cisco;  Green  River, 
Utah,  and  Bluff,  the  waters  of  the  Colorado,  Green, 
and  San  Juan  Rivers  combine  to  form  Lake  Powell. 
Many  of  the  same  changes  identified  at  these 
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Table  13.  —  Significant  differences  in  mean  concentration  and  mass  fraction  (P  <  0.05)*  between  the 
pre- 1963  and  the  post-1 965  periods,  December-March 


Station 


Calcium         Magnesium         Sodium 


Chloride 


Sulfate 


Bicarbonate 


Upper  Colorado  Basin 


Gunnison  R.  near 
Grand  Junction 

Colorado  R.  near 
Cisco 


Green  R.  near 

Greendale 
Yampa  R.  near 

Maybell 
White  R.  near 

Watson 
Green  R.  at  Green 

River,  Utah 


San  Juan  R. 
near  Bluff 


Colorado  R.  at 

Lees  Ferry 
Colorado  R.  below 

Hoover  Dam 


-59 
1 

-33 
1 


1 
-5 

-6 
(+) 


-31 


-38 

1 

-7 

-1 


-31 

-66 

'(-) 

-2 

-19 

-66 

<-) 

-1 

Green 

Basin 

- 

-19 

6 

_ 

1 

- 

- 

-16 

1 

-1 

-A 

-17 

-1 

San  Juan  Basin 

-6 

-25 

(+) 

- 

Lower  Colorado  Basin 

-17 

-65 

- 

-1 

3 

12 

(+) 

1 

-9 

-41 
2 


-5 

-2 
-22 

-3 
-17 

-2 


-59 

-2 

-18 

2 


-310 
-6 

-199 
-6 


-55 

35 
6 


-26 


-117 
-3 


-165 


-52 
7 

-31 
4 


2 

-4 

4 

-11 

2 


-30 
3 


-62 
4 


*  Upper  value  —  concentration  (mg/L),  lower  value  —  mass  fraction  (percent). 
1  (+)  or  (-)  indicate  a  change  of  less  than  0.5  percent. 


three  stations  were  also  found  at  Lees  Ferry, 
downstream  of  Glen  Canyon  Dam.  However,  the 
magnitude  of  these  changes  was  generally  greater 
than  that  which  would  have  resulted  from  a  simple 
combination  of  the  changes  upstream.  Further 
downstream,  at  the  Hoover  Dam  station,  none  of 
the  impacts  of  upstream  development  can  be 
identified.  The  two  reservoirs,  Lakes  Powell  and 
Mead,  seem  to  have  the  major  controlling  influence 
on  downstream  water  quality. 

Impacts  upstream  of  Lake  Powell  are  apparently 
minor  compared  with  the  impact  of  the  lake  itself. 
This  is  clarified  by  comparing  significant  changes 
in  the  streamf  low  and  solute  loads  throughout  the 
basin  (fig.  1 2).  Changes  in  the  load  downstream  of 
Lake  Powell  are  generally  an  order  of  magnitude 


greater  than  those  upstream.  However,  interpre- 
tation of  these  results  is  complicated  by  the 
concomitant  large  changes  in  streamflow.  Storage 
in  Lake  Powell  increased  through  most  of  the 
post-1 965  period  (the  reservoir  was  not  filled  until 
1 980).  Much  of  the  reduction  in  solute  load  may  be 
explained  by  the  overall  reduction  in  flow. 

Equally  difficult  to  explain  are  the  increases  in 
magnesium,  sodium,  and  chloride  concentrations 
identified  at  the  station  below  Hoover  Dam.  These 
increases  could  be  caused  by  the  greater  relative 
contribution  of  the  Little  Colorado  River  to  the  Lake 
Mead  inflow  during  the  filling  of  Lake  Powell.  Blue 
Springs,  near  the  mouth  of  the  Little  Colorado,  is  a 
major  source  of  sodium  chloride  loading  and  also 
contributes  magnesium,  calcium,  and  bicarbonate. 
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Table  1 4.  —  Significant  differences  in  mean  concentration  and  mass  fraction  (P  <  0.05)*  between  the 

pre-1963  and  the  post-1 965  periods,  May-June 


Station 

Calcium 

Magnesium         Sodium 

Chloride 

Sulfate 

Bicarbonate 

Upper  Colorado  Basin 

Gunnison  R.  near 

- 

- 

- 

- 

- 

Grand  Junction 

-1 

1 

- 

9 

-8 

Colorado  R.  near 

- 

5 

17 

- 

- 

Cisco 

Green  Basin 

2 

Green  R.  near 

19 

12                    29 

7 

127 

22 

Greendale 

-2 

1                       2 

- 

12 

-13 

Yampa  R.  near 

- 

3 

- 

15 

- 

Maybell 

-2 

1 

- 

4 

-6 

White  R.  near 

- 

- 

- 

- 

- 

Watson 

- 

1 

-1 

3 

- 

Green  R.  at  Green 

- 

4                        7 

- 

36 

- 

River,  Utah 

-2 

1                         1 
San  Juan  Basin 

7 

-6 

San  Juan  R. 

11 

7                     15 

5 

77 

17 

near  Bluff 

-2 

1                        1 
Lower  Colorado  Basin 

1 

8 

-8 

Colorado  R.  at 

54 

13                    49 

38 

155 

_ 

Lees  Ferry 

-3 

1     (+)                      4 

4 

11 

-16 

Colorado  R.  below 

-8 

3                       8 

16 

- 

-6 

Hoover  Dam 

-1 

(+)                       1 

2 

- 

- 

*  Upper  value  —  concentration  (mg/L),  lower  value  —  mass  fraction  (percent). 
1  (+)  or  (-)  indicate  a  change  of  less  than  0.5  percent. 


Increases  in  the  concentration  of  the  calcium  and 
bicarbonate  may  have  been  offset  by  calcium 
carbonate  precipitation  in  Lake  Mead.  The  overall 
result  would  be  the  increased  sodium  and  chloride 
concentrations  and  the  somewhat  smaller  increase 
in  magnesium,  shown  in  tables  13  through  15. 


Trends  Following  Regulation 

Significant  trends  identified  from  linear  regressions 
of  the  seasonal  concentration,  load,  and  mass 
fraction  over  time  for  the  post-1 965  period  are 
listed  in  tables  17  through  19.  Significant  trends 
are  far  less  common  than  significant  changes 
following  regulation.  This  suggests  that  the  major 
impact  of  reservoirs  on  water  quality  occurs  im- 
mediately. Where  trends  are  identified  for  regu- 


lated stations,  they  are  toward  decreasing  sea- 
sonal variability.  This  can  be  seen  by  comparing 
the  trends  in  sulfate  concentration  and  mass 
fraction  from  tables  17  through  19  with  the 
post-1 965  means  shown  in  table  1 6.  For  example, 
the  downward  trends  in  concentration  and  mass 
fraction  at  Lees  Ferry  in  May-June  will,  if  they 
continue,  reduce  the  range  of  seasonal  averages 
in  both  cases.  This  correspondence  also  occurs  for 
the  Cisco  station  (concentration  in  December- 
March,  mass  fraction  in  December-March  and 
August-November),  the  Green  River  Station  (con- 
centration and  mass  fraction  in  August-November), 
and  the  Hoover  Station  (concentration  in  December- 
March). 

Decreasing  trends  in  the  concentration  and  load  of 
sodium  and  chloride  on  the  White  River  and  the 
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Table  1 5.  —  Significant  differences  in  mean  concentration  and  mass  fraction  (P  <  0.05)*  between  the 
pre- 1963  and  the  post- 1965  periods,  August-November 


Station 

Calcium 

Magnesium 

Sodium 

Chloride 

Sulfate 

Bicarbonate 

Upper  Colorado  Basin 

Gunnison  R.  near 

-54 

-26 

-66 

-8 

-295 

_ 

Grand  Junction 

1 

- 

-2 

- 

-3 

4 

Colorado  R.  near 

-24 

-17 

-41 

- 

-188 

- 

Cisco 

1 

1     (-) 

Green 

Basin 

3 

-6 

3 

Green  R.  near 

_ 

_ 

_ 

_ 

_ 

_ 

Greendale 

- 

(+) 

-1 

- 

- 

- 

Yampa  R.  near 

- 

2 

- 

- 

21 

- 

Maybell 

- 

(+) 

- 

-2 

4 

- 

White  R.  near 

-9 

- 

- 

-22 

- 

- 

Watson 

- 

1 

- 

-3 

- 

3 

Green  R.  at  Green 

- 

- 

-7 

-11 

- 

- 

River,  Utah 

- 

(+) 

-1 

-1 

2 

- 

San  Juan  Basin 

San  Juan  R. 

-18 

_ 

_ 

_ 

-67 

_ 

near  Bluff 

- 

0.3 

- 

- 

- 

2.8 

Lower  Colorado  Basin 

Colorado  R.  at 

-29 

-18 

-61 

-39 

-223 

-55 

Lees  Ferry 

- 

- 

- 

- 

-5 

6 

Colorado  R.  below 

-7 

2 

9 

15 

- 

- 

Hoover  Dam 

-1 

(+) 

1 

2 

- 

- 

*  Upper  value  —  concentration  (mg/L),  lower  value  —  mass  fraction  (percent). 
'  (+)  or  (-)  indicate  a  change  of  less  than  0.5  percent. 


concentration,  load,  and  mass  fraction  of  sulfate  at 
Cisco  may  be  associated  with  salinity  control 
activities,  which  were  discussed  in  the  previous 
section. 

Most  of  the  increasing  trends  were  for  loads  at 
Lees  Ferry  during  the  May-June  season.  These 
were  associated  with  an  increasing  trend  in 
streamflow,  the  only  flow  trend  identified.  Because 
no  solute  concentration  trends  were  found,  the 
flow  trend  is  probably  the  primary  reason  for  the 
increasing  loads.  Overall,  the  impacts  of  reservoir 
storage  are  easily  identifiable  using  the  two 
selected  statistical  tests.  Other  activities  influ- 
encing water  quality  in  the  basin  were  also  found. 
However,  the  reservoirs,  especially  Lakes  Powell 
and  Mead,  appear  to  exert  major  control  on  the 
quality  of  the  inflow  to  the  lower  basin. 


CONCLUSIONS 

1 .  A  procedure  was  proposed  and  FORTRAN  pro- 
grams were  developed  for  computing  monthly 
values  for  streamflow,  concentrations  of  six  major 
ions,  and  concentration  of  the  sum  of  these  six 
ions.  This  procedure  was  based  on  USGS 
WATSTORE  data  of  daily  flows,  daily  specific 
conductance,  and  approximately  monthly  chemical 
analyses.  Monthly  values  were  generated  for  15 
gaging  stations  in  the  Colorado  River  Basin  after 
deletion  of  obviously  erroneous  WATSTORE  data 
and  identification  of  suspected  erroneous  data. 
The  resultant  computer  files  provide  an  easily 
accessible  historical  record  of  flow  and  salinity, 
which  can  be  used  in  statistical  analyses  by 
Government,  academic,  and  private  researchers. 
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Table  1 7.  —  Significant  trends  in  the  post-1 965  period,  December-March 


Linear  trend  slope,  a'  (P  <  0.05)* 


Station 


Flow 
m3/s 


Calcium     Magnesium     Sodium        Chloride        Sulfate     Bicarbonate 


Gunnison  R.  near 
Grand  Junction 


Colorado  R.  near 
Cisco 


Green  R.  near 
Greendale 


Yampa  R.  near 
Maybell 


Upper  Colorado  Basin 


0.07 


- 

-2.66 

-0.02 

-0.08 

-0.64 

_ 

-0.05 

-0.13 

-0.03 

- 

Green  Basin 

-0.51 


-5.77 
-0.40 
-0.25 


-0.35 
-0.07 


0.22 


0.79 


White  R.  near 
Watson 


-0.04 


-3.34 
-0.06 


-3.68 
-0.06 


-0.09 


Green  R.  at  Green 
River,  Utah 


San  Juan  R. 
near  Bluff 


-0.07 


-0.54 
-0.06 
-0.07 


San  Juan  Basin 


1.54 
0.25 


Colorado  R.  at 
Lees  Ferry 


Colorado  R.  below 
Hoover  Dam 


Lower  Colorado  Basin 
-0.89 

-0.04  0.02 

-0.61  -  -0.82 


-1.16 


-2.25 
-0.07 


*  Upper  value  —  concentration  (mg/L) 
middle  value  —  load  (kg/s) 
lower  value  —  mass  fraction  (percent). 

After  the  proposed  procedure  is  reviewed,  modified, 
and  approved,  the  data  f  iies  will  provide  a  common 
basis  for  the  evaluation  of  salinity  relationships  in 
the  Colorado  River  Basin.  Therefore,  discussions 


of  differing  statistical  results  regarding  salinity 
relationships  will  center  on  differences  in  the 
analyses  themselves  and  not  on  differences  in  the 
data. 
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STREAMFLOW  (mVt) 


CHLORIDE    (kg/*) 


o  =  gages  station  location 
(see  fig.  1  for  identification) 


SULFATE  (kg/*) 


BICARBONATE  (kg/*) 


Figure  14.  —  Significant  differences  (P<  0.05)  in  streamflow(mVs)  and  anion  loads  (kg/s)  between  the  pre- 1963 
and  the  post-1965  periods.  Upper  value:  December-March,  middle  value:  May-June,  lower  value:  August- 
-November. 
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Table  1 8.  —  Significant  trends  in  the  post-1965  period,  May-June 


Linear  trend  slope,  a"1  (P  <  0.05)' 


Station  Flow  Calcium     Magnesium     Sodium        Chloride        Sulfate     Bicarbonate 

m3/s 


Upper  Colorado  Basin 

Gunnison  R.  near  _______ 

Grand  Junction  ______ 

Colorado  R.  near  _______ 

Cisco  ______ 

-0.07  _  _  _  _ 

Green  Basin 

Green  R.  near  -  -  -0  39  -  -0.41 

Greendale  ______ 

-0.03  -  -0  05  -  0.17 

Yampa  R.  near  _______ 

Maybell  ______ 

0.40 

White  R.  near  _______ 

Watson  -  -0.05 

*-0.30 

Green  R.  at  Green  _______ 

River,  Utah  _  _  _  -0.14 

San  Juan  Basin 

San  Juan  R.  _______ 

near  Bluff  -  -0.14 

Lower  Colorado  Basin 

Colorado  R.  at  -1.07  -  -  -  -2.90 

Lees  Ferry  -0.91  -  -2.81 

-0.06  -  -0.10 

Colorado  R.  below  -  -0.48  -0.19  -0.80  - 

Hoover  Dam  ______ 


*  Upper  value  —  concentration  (mg/L) 
middle  value  —  load  (kg/s) 
lower  value  —  mass  fraction  (percent). 

Based  on  the  monthly  values  of  concentrations  2.  A  theoretical  residence  time  model  for  describ- 

and  streamflow,  the  annua!  cycle  appears  to  have  ing  the  concentration-streamflow  relationship  at 

three  distinct  seasons:  baseflow  (December-  each  station  with  unregulated  flows  was  derived 

March),  peak  runoff  (May-June),  and  return  flow  and  compared  with  the  widely  used  power  model. 

(August-November).  For  monthly  values  of  concentration  and  stream- 
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Table  1 9.  —  Significant  trends  in  the  post-1965  period,  August-November 


Linear  trend  slope,  a'(P^  0.05)* 


Station  Flow  Calcium     Magnesium     Sodium        Chloride        Sulfate     Bicarbonate 

m3/s 


Upper  Colorado  Basin 


Gunnison  R.  near 
Grand  Junction 


0.09 


-0.06 
-0.02 


-0.10 


-0.07 


-0.44 


Colorado  R.  near 
Cisco 


-0.13 


-0.53 
-0.27 


Green  R.  near 
Greendale 


Upper  Green  Basin 
-0.32 


-0.26 


Yampa  R.  near 
Maybell 


-0.01  -0.01  -0.01 


White  R.  near 
Watson 


-0.03  -0.03 


Green  R.  at  Green 
River,  Utah 


-0.66  -0.33  -0.70 


-3.09 

-0.15  -0.20 


San  Juan  Basin 


San  Juan  R. 
near  Bluff 


-2.33 


-0.15 


-0.52 


-3.42 


Lower  Colorado  Basin 


Colorado  R.  at 
Lees  Ferry 


9.85 


0.59  0.26  0.73  0.55  2.28  1.41 

-0.04  - 


Colorado  R.  below 
Hoover  Dam 


-0.56 


-0.89 


*  Upper  value  —  concentration  (mg/L) 
middle  value  —  load  (kg/s) 
lower  value  —  mass  fraction  (percent). 
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flow,  the  models  describe  the  relationship  equally 
well,  and  both  can  be  used  to  correct  for  the 
variability  of  streamflow  in  evaluating  the  statis- 
tical significance  of  long-term  changes  in  con- 
centration. 


load  increased  by  0.4  percent  per  year  (5  * 
106  kg/a),  and  the  bicarbonate  load  increased  by 
0.7  percent  per  year  (21  *  106  kg/a).  No  statis- 
tically significant  changes  in  the  chloride  load  are 
evident. 


The  relationships  between  the  solute  concentra- 
tions and  the  streamflow  vary  with  the  solute, 
season,  and  gaging  station.  The  bicarbonate  ion 
concentration  is  relatively  independent  of  the 
streamflow  for  all  stations.  Because  the  concen- 
tration-streamflow  relationship  varies  from  solute 
to  solute,  TDS-streamflow  equations  may  be  less 
accurate  than  equations  for  individual  solutes. 

3.  Long-term  trends  in  solute  concentrations  for 
12  gaging  stations  were  evaluated  for  statistical 
significance  by  three  empirical  equations  with 
multiple  linear  regressions  and  by  the  residence 
time  equation  with  a  nonlinear  regression.  These 
four  equations  yielded  similar  results  and  are, 
therefore,  corroborative. 

Long-term  trends  in  solute  loads  can  be  estimated 
from  the  regression  results  for  concentrations. 
With  a  linear  streamflow  term  in  the  regression, 
this  procedure  is  preferable  to  direct  regressions 
with  solute  loads  because  upstream  diversions 
have  a  smaller  effect  on  the  estimated  trends. 
With  a  power  streamflow  term,  both  procedures 
yield  identical  results,  making  load  regressions 
redundant. 


4.  After  construction  of  the  Colorado  River  Storage 
Project  reservoirs  in  the  early  1 960's,  the  variability 
of  solute  concentrations  and  mass  fractions  were 
reduced  at  many  downstream  gaging  stations.  The 
degree  of  reduction  was  relatively  proportional  to 
the  degree  of  regulation.  At  some  sites,  trends 
toward  further  reduction  in  variability  developed 
after  the  initial  reservoir  filling.  Several  other 
activities  of  man,  including  mining,  urbanization, 
and  salinity  control  practices,  may  have  affected 
the  water  quality  in  the  basin.  Impacts  decreased 
with  distance  downstream,  but  were  often  identi- 
fied at  more  than  one  gaging  station.  However,  the 
major  reservoirs,  Lakes  Powell  and  Mead,  masked 
all  impacts  occurring  in  the  Upper  Basin. 

No  significant  trends  in  the  load  for  any  solute 
have  been  identified  for  the  station  below  Hoover 
Dam  since  1 965;  however,  the  concentration  of  all 
solutes  except  chloride  decreased  in  at  least  one 
season. 


IMPLICATIONS  FOR  THE  SALINITY 
CONTROL  PROGRAM 


Long-term  trends  in  solute  mass  fractions  for  12 
gaging  stations  were  evaluated  for  statistical  sig- 
nificance by  three  empirical  equations  with  multiple 
linear  regression.  As  expected,  mass  fraction  is 
the  least  sensitive  to  variability  of  streamflow,  so 
that  both  the  linear  and  power  corrections  for 
streamflow  were  satisfactory.  The  results  from  the 
three  equations  were  almost  identical.  These 
results  agreed  with  the  estimated  trends  in  con- 
centration, but  yielded  no  new  information. 

Before  construction  of  the  major  upstream  dams, 
the  concentrations  of  magnesium,  sodium,  and 
sulfate  ions  during  the  basef  low  season  of  Decem- 
ber-March decreased  in  stoichiometric  propor- 
tions, resulting  in  a  reduction  of  approximately 
2  mg/L  per  year  in  the  baseflow  TDS  of  the 
Colorado  River  at  Lees  Ferry  and  at  the  Grand 
Canyon.  At  the  Grand  Canyon  from  1 925-61 :  the 
annual  magnesium  load  decreased  at  an  average 
rate  of  0.4  percent  per  year  (2  *  106  kg/a),  the 
sodium  and  sulfate  loads  decreased  0.3  percent 
per  year  (5  *  1 06  and  1 2  *  1 06  kg/a),  the  calcium 


The  development  of  a  standard  data  base  for 
salinity  at  the  major  gaging  stations  in  the  Colorado 
River  Basin  will  benefit  researchers  evaluating  the 
salinity  problem  and  potential  measures  for  salinity 
control.  The  present  work  indicates  that  standard 
statistical  procedures  applied  to  these  data  are 
sensitive  enough  to  identify  impacts  of  the  mag- 
nitude projected  for  proposed  USBR  salinity  control 
projects.  However,  it  also  shows  the  overriding 
influence  of  Lakes  Powell  and  Mead  on  salinity  in 
the  Colorado  River  downstream  of  Hoover  Dam. 
More  research  is  needed  to  understand  the 
processes  occurring  within  these  reservoirs  to 
assess  the  probable  impact  of  the  Upper  Basin 
projects  on  Lower  Basin  water  quality. 

Salinity  downstream  of  Lake  Mead  was  found  to 
be  relatively  stable,  with  some  recent  decreases  in 
concentration.  Therefore,  the  observed  changes 
further  downstream  at  Imperial  Dam  may  have 
been  caused  by  changing  conditions  in  the  Lower 
Basin  area.  Additional  studies  are  needed  to 
identify  their  sources. 
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RECOMMENDATIONS  FOR  FURTHER 
RESEARCH 

After  a  review  and  any  necessary  revisions,  the 
data  base  developmenttechnique  should  be  applied 
to  other  major  stations  in  the  basin,  particularly 
stations  on  the  Colorado  River  below  Parker  Dam 
and  at  Imperial  Dam.  In  addition,  the  existing  data 
base  should  be  updated  as  new  data  become 
available.  Potential  modifications  to  the  technique 
include  increasing  the  period  of  interpolation  for 
chemical  analyses,  particularly  for  stations  with 
regulated  flow,  where  the  chemical  variability  is 
small.  This  would  extend  the  period  of  record  at 
sites  where  the  USGS  has  gone  to  2-month  or 
quarter-year  sampling  intervals. 

Additional  statistical  analyses  are  suggested  using 
this  expanded  data  set.  Seasonal  division  of  the 
data  for  stations  downstream  of  major  impound- 
ments is  not  necessary  because  of  their  reduced 
variability.  Regression  analyses  should  be  con- 
ducted only  for  concentration,  with  meq/L  as  the 
preferred  units.  This  would  facilitate  the  identi- 
fication of  chemical  relationships  in  significant 
changes  and  trends.  Differences  in  the  load 
between  consecutive,  unregulated  stations  should 
also  be  analyzed  to  identify  impacts  such  as  those 
hypothesized  for  the  Grand  Valley  area. 

The  greatest  need  is  for  a  more  complete  investi- 
gation of  the  potential  sources  of  identified  changes 
and  trends.  Statistical  results  should  receive  a 
multi-disciplinary  review  to  determine  sources.  If 
possible,  data  should  be  collected  to  confirm 
suspected  sources.  The  probable  impact  of  USBR 
salinity  control  projects  can  then  be  evaluated. 
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GENERATION  OF  MONTHLY 

AVERAGE  STREAMFLOW  AND 

CONCENTRATION  VALUES 

A  series  of  FORTRAN  computer  programs  was 
developed  to  check  data  retrieved  from  the  USGS 
WATSTORE  data  base  and  then  to  use  these  data 
to  generate  average  monthly  values  for  discharge 
and  constituent  concentrations.  These  programs 
are  described  below  in  the  order  of  their  application 
to  a  particular  set  of  retrieved  data.  Also  described 
is  the  iterative  procedure  used  to  identify  possible 
errors  in  the  WATSTORE  data.  The  programs  are 
stored  on  the  USBR  CYBER  under  User  No. 
ER1 000D.  Catalog  names  are  identical  to  program 
names. 

Modifications  made  to  the  WATSTORE  data  before 
generating  monthly  values  are  tabulated  for  the 
appropriate  station  in  table  2  of  appendixes  B 
through  P.  The  responsible  USGS  district  offices 
are  notified  of  discrepancies  between  WATSTORE 
data  and  published  Water  Supply  Paper  values. 


Step  1  —  WATSTORE  Data  Retrieval 

Data  were  transferred  from  the  USGS  computer 
system  in  Reston,  Virginia,  to  the  USBR  CYBER  by 
an  automatic  communication  link  developed  by 
Main  [1  ]*.  Three  retrievals  were  required  for  each 
station:  (1 )  daily  value  streamflow,  (2)  daily  value 
specific  conductance,  and  (3)  periodic  water 
quality.  Specific  conductance  retrievals  included 
both  random  daily  measurements  (statistic  code 
1 1)  and  mean  daily  values  (statistic  code  3).  The 
latter  are  available  only  for  stations  with  continuous 
recorders.  Water  quality  parameters  retrieved  were 
calcium  (code  91 5),  magnesium  (code  925),  sodium 
(codes  930),  sodium  plus  potassium  (code  933), 
chloride  (code  940),  sulfate  (code  945),  bicarbonate 
(code  440,  90440,  95440,  and  99440),  alkalinity 
(codes  410,  90410),  specific  conductance  (code 
95),  andTDS  (total  dissolved  solids)  as  evaporation 
residues  at  180  °F  (code  70300).  In  all  cases  the 
entire  period  of  record  was  retrieved. 


Step  2  —  Water  Quality  Inventory 

The  water  quality  data  file  can  be  inventoried  using 
program  QWTAB1.  Output  from  this  program 
gives  the  period  of  record,  number  of  samples,  and 


*  Numbers  in  brackets  refer  to  entries  in  the  bibliography. 


range  of  values  for  each  parameter  in  the  file  (fig. 
A-1).  Values  which  appear  too  extreme  can  be 
located  in  the  data  file,  checked  against  the 
appropriate  water  supply  paper,  and  corrected,  if 
necessary.  In  this  and  all  subsequent  programs,  a 
zero  water  quality  value  is  considered  missing 
data. 

The  inventory  is  generated  from  an  interactive 
terminal  using  the  following  commands: 

BAT 

GET,  QWTAB1,  data  file 
FTN,  l  =  QWTAB1,L  =  0 
LGO,  data  file 

where  "data  file"  is  the  catalog  name  of  a 
WATSTORE  water  quality  data  file.  Output  will  be 
displayed  on  the  terminal.  A  separate  output  file 
on  TAPE4  lists  duplicate  and  overlapping  dates. 
This  list  can  be  displayed  on  the  terminal  by 
entering: 

REWIND,  TAPE4 
LIST,F=TAPE4 


Step  3  —  TDS  and  Conductivity  Regressions 

Three  programs  (REGWQ1,  REGWQ2,  and 
REGWQ3)  are  used  to  regress  the  sum  of  six 
solutes  against  the  evaporation  residue  TDS  and 
to  determine  the  empirical  constant  in  the  relation- 
ship between  the  observed  and  computed  specific 
conductance. 

The  first  program  (REGWQ1)  reads  the  WATSTORE 
data  and  computes  the  dependent  and  indepen- 
dent regression  variables  for  all  samples  with 
sufficient  data.  If  sodium  concentration  (code  930) 
is  missing,  sodium  plus  potassium  data  (code  933) 
are  used.  If  bicarbonate  concentration  is  missing, 
it  is  computed  from  alkalinity: 

Bicarbonate  =  alkalinity /0.8202 

This  is  simply  the  mass  ratio  between  bicarbonate 
(HCO3)  ar,d  calcium  carbonate  (CaC03)  in  the 
stoichiometric  relationship: 

CaC03(s)  +  H20  +  C02(g)  «*  2HC031  +  Ca+2 

The  conversion  is  made  under  the  assumption  that 
the  total  alkalinity  (as  CaC03)  results  from  bicar- 
bonate. This  assumption  is  reasonable  for  water 
with  pH  less  than  8.3  [2].  Because  bicarbonate  and 
alkalinity  analyses  in  WATSTORE  are  specified 
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83/08/03.    UATER  QUALITY  DATA  FILE  INVENTORY 

FOR 
STATION  -  9180500 
COLORADO  RIVER  NEAR  CISCO  UTAH 


15.42.40. 


NO.  OF 

PARflETI 

FIRST 

LAST 

NO.  OF 

FLAGD 

ID 

CODE 

DATE 

DATE 

ITEMS 

minimum 

MAXIMUM 

MEAN 

ITEMS 

1. 

ALK(F) 

410 

2810 

8009 

1237 

80.00 

221.00 

152.97 

0 

2. 

HC03CU) 

440 

2810 

8006 

1345 

97.00 

270.00 

185.43 

0 

3. 

CA 

915 

2810 

8305 

1304 

36.00 

265.00 

111.57 

0 

4. 

flG 

925 

2810 

8305 

1303 

6.30 

131.00 

44.70 

0 

5. 

NA 

930 

2810 

8305 

1154 

15.00 

495.00 

158.22 

0 

6. 

CL 

940 

2810 

8305 

1374 

8.00 

695.00 

139.43 

0 

7. 

S04 

945 

2810 

8305 

1371 

52.00 

1320.00 

430.36 

0 

8. 

SC 

95 

3510 

8305 

1363 

330.00 

4049.99 

1442.13 

0 

9. 

TDS(R) 

70300 

3510 

8304 

976 

208.00 

2350.00 

972.13 

0 

10. 

NA+K 

933 

4310 

8002 

208 

21.00 

298.00 

127.32 

0 

11. 

ALK(L) 

90410 

8010 

8305 

24 

83.00 

180.00 

141.96 

0 

Figure  A-1 .  —  Example  output  from  program  QWTAB1 


from  several  sources,  the  determination  is  at- 
tempted in  the  following  order:  (1)  bicarbonate, 
titration  to  pH  4.5,  lab  (parameter  code  95440),  (2) 
bicarbonate,  unspecified  (code  440),  (3)  bicar- 
bonate, incremental  titration,  lab  (code  90440),  (4) 
alkalinity,  titration  to  pH  4.5,  lab  (code  9041 0),  (5) 
bicarbonate,  incremental  titration,  field  (code 
99440),  and  (6)  alkalinity,  titration  to  pH  4.5,  field 
(code  41 0).  In  the  study  area,  only  sources  (2),  (4), 
and  (6)  are  generally  available. 

To  compare  the  sum  of  six  solutes  with  evaporation 
residue  TDS,  the  residue  TDS  must  be  increased  to 
account  for  bicarbonate  lost  during  evaporation  at 
180  °F: 


2HC03  +  heat  =*  C032  +  H20(g)  +  CO 


'2(g) 


According  to  the  mass  ratio  of  C032  to  2HC03, 
50.83  percent  of  the  bicarbonate  is  lost.  Therefore, 
the  regression  equation  is 

sum  =  a  +  b  {TDS(res,due)  +  0.5083  (HC03)} 

The  sum  of  six  solutes  (sum)  and  the  adjusted  TDS 
are  computed  by  program  REGWQ1  for  each 
complete  analysis.  Regression  parameters  a  and  b 
are  evaluated  using  program  REGWQ2.  Expected 
values  of  TDS  regression  coefficients  are  a  slope 
(b)  of  approximately  1  and  a  slightly  negative 
intercept  (a).  This  accounts  for  constituents,  mainly 
silica  and  potassium,  included  in  the  evaporation 
residue  but  not  in  the  sum  of  major  ions. 


Specific  conductance  (L)  can  be  estimated  from 
the  equivalent  ionic  conductances  using  the 
method  of  Tanji  and  Biggar  [3]: 

L  =  %(k°i-aCp5)Ci 


where:     A°,    =  Limiting  ionic  conductance  of 
species  i 
C,     =  Molar  concentration  of  species  i 
a      =  An  empirical  constant 

Equation  5  may  be  rearranged: 

2,(\0iC,)-L  =  aC>s 


or: 


y-  ax 


REGWQ1  computes  the  variable  (x  and  y)  terms  for 
each  complete  analysis  in  the  WATSTORE  water 
quality  data  file.  The  constant  a  is  then  evaluated 
using  linear  regression  forced  through  the  origin 
(program  REGWQ2).  The  normal  range  of  values 
for  a  is  5-1 0. 

Linear  regressions  are  computed  using  SPSS 
procedures  [4].  Program  REGWQ2  contains  the 
SPSS  directives  to  produce  line-printer  scatter 
diagrams  of  the  data  (output  from  REGWQ1 )  and 
the  linear  regression  coefficients.  Examples  are 
shown  in  figures  A-2  and  A-3. 
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Figure  A-2.  —  Example  output  from  regression  of  sum  of  solutes  on  adjusted  residue  TDS. 
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Figure  A-3.  —  Example  output  from  regression  evaluation  of  conductance  constant  a. 
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By  inspection  of  the  scatter  diagrams,  obvious  data 
problems  can  be  identified.  These  are  located  in 
the  data  file,  checked  against  Water  Supply  Papers 
and  corrected.  Final  regression  coefficients  are 
used  in  later  programs. 

The  third  program  (REGWQ3)  is  the  Batch  JCL  (job 
control  language),  which  runs  the  other  two. 
Before  submission,  it  must  be  edited  to  access  the 
desired  water  quality  input  file  and  create  the 
desired  output  file.  The  program  is  then  run,  using 
the  commands: 

OLD,  REGWQ3 
SUB 


conductance  and  streamflow,  outside  which  ex- 
treme values  and  dates  will  be  listed. 

The  REGSC3  output  is  also  used  to  find  the  slope  of 
the  log-log  relationship  between  L  and  Q.  This 
slope  is  identical  to  the  exponent,  B,  in  the 
relationship: 

L=aQB 

The  fitted  exponent  is  used  in  step  6  to  determine 
the  weighting  factor  for  interpolation  of  missing 
conductance  data. 


A  dayfile,  which  includes  cost  information,  will  be 
saved  under  the  name  DFRTDS. 


Step  4  —  Regression  of  Conductivity  Versus 
Streamflow 

The  daily  value  data  for  specific  conductance  (L) 
and  streamflow  (Q)  are  also  checked  by  regres- 
sion. As  in  step  3,  a  series  of  programs  is  used. 
REGSC1  reads  data  from  the  two  files,  sorts  them 
by  date,  and  transforms  them  to  natural  logarithms. 
Only  dates  with  both  an  L  and  Q  value  are  used  in 
the  regression.  REGSC2  is  the  SPSS  directive  file. 
It  produces  a  scatter  diagram  and  regression 
statistics  (fig.  A-4).  REGSC3  is  the  CYBER  job 
control  file  that  executes  the  other  two. 

To  run  the  procedure,  REGSC3  must  be  edited  to 
access  the  desired  data  files.  An  output  file  name 
can  also  be  specified.  From  an  interactive  terminal, 
enter  the  following  commands: 

OLD,  REGSC3 

(edit) 

SUB 

Output  will  be  sorted  in  the  specified  file.  A  dayfile 
will  also  be  created  under  the  name  DFRSC. 

The  scatter  diagram  can  be  used  to  identify  out- 
liers, which  may  be  data  entry  errors.  The  dates  of 
outliers  and  extreme  values  can  be  listed  using 
program  CHKSC.  CHKSC  is  run  in  BATCH  mode 
using  the  FTN  and  LGO  commands.  Before  it  is 
run,  the  program  must  be  edited  to  specify  the 
standard  error  of  estimate  from  REGSC3  output. 
Also,  the  user  must  set  the  ranges  of  specific 


Step  5  —  Tabulation  of  Water  Quality  Data 

The  WATSTORE  water  quality  data  file  is  re- 
formatted in  a  readable  table  by  program  QWTAB2. 
If  necessary,  bicarbonate  and  sodium  values  are 
determined  as  described  in  step  2.  Single  missing 
concentrations  within  any  individual  sample  are 
then  estimated  assuming  a  charge  balance  of  zero. 
Estimates  are  constrained  to  values  within  the 
observed  range  for  that  constituent.  If  two  or  more 
values  are  missing  from  a  single  analysis,  no 
estimation  is  attempted. 

QWTAB2  also  checks  the  estimates  of  sum  of 
constituents  from  residue  TDS,  the  calculation  of 
specific  conductance  from  individual  ionic  equi- 
valent values,  and  the  charge  balance.  The  TDS 
and  conductance  estimates  are  made  using  the 
regression  coefficients  evaluated  in  step  3.  The 
charge  balance  is  computed  using  the  following 
factors  to  convert  concentrations  from  mg/L  to 
meq/L: 


Constituent 

Factor 

Calcium 

20.04 

Magnesium 

12.15 

Sodium 

22.99 

Chloride 

35.45 

Sulfate 

48.03 

Bicarbonate 

61.02 

A  partial  listing  of  the  output  table  is  shown  in 
figure  A-5. 

A  separate  table  (fig.  A-6)  gives  error  statistics  for 
the  computed  TDS,  specific  conductance,  and  net 
charge. 
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Figure  A-4.  —  Example  output  from  log-log  regression  of  daily  value  conductance  on  streamflow. 
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Deviations  between  observed  and  estimated  values 
are  accumulated  and  used  to  compute  two  error 
statistics,  RMSE  (root-mean-square  error)  and 
NME  (normalized  mean  error): 


values  generation  program  in  step  6.  After  final 
data  corrections,  the  WATSTORE  water  quality 
data  file  is  no  longer  needed  and  may  be  purged  or 
archived. 


NME 


RMSE  = 


100 


=  [ 


X  (X-X)2 


10.5 


vr  abs(x-x)  i 


where:     X    =  Observed  value  (zero  for  charge 
balance) 
X    =  Estimated  value 
n    =  Number  of  samples 

RMSE  is  in  the  units  of  X;  NME  is  a  percentage 
value. 

To  facilitate  error  analysis,  the  dates  with  signifi- 
cant deviations  between  measured  and  predicted 
sum  of  constitutents,  specific  conductance,  and 
charge  balance  (assuming  the  normal  charge 
balance  is  zero)  are  also  listed  with  the  error  table. 
Data  for  these  dates  can  be  checked  by  comparison 
with  USGS  Water  Supply  Papers.  If  specific  errors 
are  not  found,  unreasonable  values  may  be  set  to 
zero  and  new  values  will  be  estimated  by  the 
program.  If  changes  are  made,  steps  3,  4,  and  5 
must  be  reprocessed. 

QWTAB2  requires  that  site-specific  data  be  initial- 
ized in  the  program.  These  data  include  the  TDS 
regression  coefficients  and  the  conductance  con- 
stant (a)  from  step  3,  and  the  log-log  conductance 
versus  flow  regression  slope  from  step  4.  Two 
dates  are  also  required,  the  beginning  of  specific 
conductance  monitoring  and  the  beginning  of  flow 
regulation  by  an  upstream  dam.  The  former  is 
determined  by  the  first  date  of  daily  specific 
conductance  under  statistic  code  3.  Flow  regulation 
dates  are  generally  given  under  "Remarks"  in  the 
station  description  in  Water  Supply  Papers. 

Site-specific  values  are  set  using  XEDIT.  The 
source  code  must  tnen  be  replaced  to  save  the 
edited  values.  The  program  is  run  by  submitting 
procedure  file  RQWT2.  This  must  also  be  edited  to 
read  the  desired  input  and  to  specify  output  file 
names  for  both  the  data  and  error  analysis  tables. 

The  data  and  error  files  in  disk  storage  may  be 
printed  on  the  line  printer  or  Xerox.  The  error  file 
may  then  be  purged,  but  the  data  file  must  be 
preserved,  because  it  is  input  to  the  monthly 


Step  6  —  Monthly  Value  Generation 

Program  QWMON  uses  results  from  the  previous 
five  steps  to  generate  the  mean  monthly  values  for 
streamflow,  conductance,  and  solute  concentra- 
tions. These  are  computed  from  data  in  the  three 
WATSTORE  data  files  and  previous  estimates  for 
missing  values.  The  program  is  divided  into  three 
parts:  (1 )  compilation  of  daily  data,  (2)  estimation  of 
missing  values  between  samples,  and  (3)  compu- 
tation of  monthly  means. 

Compilation  of  daily  data.  —  QWMON  first  as- 
sembles data  from  three  sources  into  a  single  daily 
values  file.  The  sources  are  the  two  WATSTORE 
daily  values  files  for  streamflow  and  specific 
conductance,  and  the  QWTAB2  output  which 
includes  WATSTORE  water  quality  data  and  esti- 
mates for  missing  values  within  samples.  This  part 
of  the  program  simply  increments  the  Julian  date 
(based  on  3  366-day  calendar)  and  assigns  to  a 
scratch  file  the  proper  daily  values  for  flow,  con- 
ductance, solute  concentrations,  and  sum  of  six 
ions.  Previously  estimated  values  are  flagged.  If  an 
analysis  is  incomplete  and  the  sum  of  six  ions 
cannot  be  computed,  the  estimate  using  residue 
TDS  will  be  used  and  flagged.  Missing  values  for 
streamflow  are  labeled  "-1 "  to  differentiate  them 
from  true  zero  flow  values.  Water  quality  analyses 
for  composite  samples  are  assigned  to  each  day 
within  the  composite  period.  Discrete  sample 
analyses  are  assigned  only  to  the  day  the  sample 
was  collected. 

Estimation  of  missing  values.  —  Following  initial 
compilation,  missing  daily  values  are  filled  in  by 
weighted  linear  interpolation,  as  described  in  the 
body  of  this  report.  The  actual  mechanics  of 
estimating  missing  values  involves  creating  a 
partial  array  of  the  data  and  computing  values 
where  zeros  or  -1  's  (for  streamflow)  occur.  Array 
element  numbers,  rather  than  Julian  dates,  define 
the  interpolation  period  to  avoid  problems  with 
interpolating  over  New  Year's  Day  or  February  29. 
The  array  is  sized  so  that  all  missing  values  in  the 
first  61  elements  will  be  estimated  if  sufficient 
data  are  available.  To  ensure  this,  81  additional 
elements  must  be  included,  making  an  array  of 
1 42  elements.  This  sizing  depends  on  the  selection 
of  maximum  interpolation  intervals,  as  defined 
below.  After  completion  of  the  computational  loop, 
the  first  60  elements  are  written  to  a  scratch  tape; 
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Figure  A-5.  —  Partial  data  table  output  from  program  QWTAB2. 
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ESTIMATED  VALUE  ERROR  TABULATION 
STATION   9180500  -  COLORADO  RIVER  NEAR  CISCO  UTAH 


FURORS 


PARAMETER 


NUMBER  OF 
OBSERVATIONS 


MF  AN 
(OI?S) 


VAI  UF.S 
(PRtD) 


MAXIMUM 
(VALUE)  (DATE) 


NORMALIZED  MEAN 
(PERCENT) 


ROOT  MEAN-SQUARE 


TDS  (MG/L)  05R        1037.7   1037.7 

SPECIFIC  CON  (US/CM)     1011        1415.3   1450  8 
NET  CHARGE  (MEO/L)      1300         .060      0.0 


235?  650914 
644  O  651 105 
1.279    7804  19 


2.0 

4.0 

.6 


28.  2 
74.9 
.245 


ALPHA 


6.  28 


BETA  1 


-  13.87 


BETA2  = 


97559 


DIFFERENCE  BETWEEN  OBSERVED  AND  PREDICTED  VALUES 
(+  =  SIGNIFICANT  AT  95%  LEVEL  OE  CONFIDENCE) 


DATE 


NET  CHARGE 
(MEO/L) 


DISSOLVED  SOLIDS 
(MG/L) 


SPECIFIC  COND 
(US/CM) 


29032  1 
30102 1 
301 101 
30  1  1  1  1 
3012  1  1 
310101 
310105 
3101 1  1 
310221 
33082  1 
3401 1 1 
34072 1 
340801 

3408  1  I 

3409  1  1 
38 101  1 
38 10?  1 
501012 
501017 
50  10? 1 
501 101 
501201 
5012 1 1 
5102  10 
510220 
51032 1 
5  10901 
51092 1 
51 102 1 

52  102  1 
521 101 
530920 

53  1001 
53  1009 


.  523 
.605 
.577 
.669 
.481 
.515 
.515 
.512 
.517 
.593 
.499 
.565 
.606 
.52  1 
.553 
.020 
.022 
.346 
.34<5 
.030 
.067 
.665 
.  703 
.531 
.708 
.  488 
.  481 
.550 
.572 
.330 
.  239 
.594 
.538 
.538 


O 
O 
0 
O 
O 
O 
O 
O 
O 
O 
O 
O 
O 
O 
O 
-58 
-61  . 
-55. 
-55. 
-56, 
-74  . 
-44 
-40 
-  19 
-7. 
-28 
32 
-41  . 
-29. 
55. 
-59 
103. 
24. 
24. 


O.O 
0.0 
0.0 
0.0 


O.O 

O.O 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
-82.8 
-60.9 

-  17. 

-  17. 

3. 

43. 
-44  . 
-44 

18. 
-20. 

22. 
-97. 
-78 
-95. 
-28. 
-2?  . 

-5. 

13. 

13. 


Figure  A-6.  —  Example  of  error  statistics  table  from  program  QWTAB2. 
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indexes  (k)  on  the  remaining  82  elements  are 
redefined  to  k-60;  and  the  next  82  daily  values  are 
read  into  the  bottom  of  the  array. 

Computation  of  monthly  means.  —  Measured  and 
estimated  daily  values  are  then  used  to  compute 
the  mean  monthly  streamflow  (m3/s)  and  solute 
concentrations  (mg/L).  Output  is  written  to  TAPE7 
and  stored  on  a  user-specified  disk  file. 

To  aid  in  judging  the  quality  of  the  generated  data, 
a  second  output  file  is  written,  which  tabulates  the 
percentage  of  days  with  estimated  values  for  each 
month  and  each  parameter.  Monthly  mean  solute 
concentrations  based  on  a  high  percentage  of 
estimates  may  be  inaccurate,  especially  when 
much  of  the  conductance  data  are  also  estimated. 
This  may  be  justification  for  the  exclusion  of  data 
points  identified  as  outliers  in  later  analyses. 

Operating  procedure.  —  QWMON  is  run  in  BATCH 
mode  by  submitting  a  procedure  file,  RQWMON. 
Prior  to  submittal,  the  following  lines  must  be 
edited  as  indicated,  to  access  the  desired  input 
files  and  to  specify  the  output  file  names. 

GET,  TAPE1  =  Output  from  QWTAB2, 
TAPE2  =  streamflow  data, 
TAPE3  =  conductance  data 
1,  CALL,  SAVER(l=TAPE6,J=estimate  table 

file  name) 
CALL,  SAVER(l=TAPE7,J=output  data  file 
name) 

In  addition  to  the  monthly  means  and  percentage 
estimates  tables,  this  procedure  produces  an 
output  file,  OQWM,  which  includes  messages 
written  during  program  execution,  run  cost  data, 
and  a  complete  dayfile. 


Step  7  —  Plotting  Monthly  Values 

Generated  monthly  data  are  plotted  to  check  for 
reasonable  relative  values.  Program  PLTWQM 
plots  streamflow  and  solute  concentration  values 
from  QWMON  output  and  uses  flow  and  concen- 
tration to  compute  and  plot  the  solute  load  (kg/s). 
Because  plots  are  made  on  10-year  intervals, 
several  plots  may  be  needed  to  display  the  histori- 
cal record  for  any  single  parameter.  The  cations 
and  anions  are  grouped  on  common  plots  to  show 
relative  values.  Monthly  data  points  are  connected 


for  sequential  months.  Possible  errors  may  be 
determined  visually  and  checked  in  intermediate 
output  tables. 

The  program  is  run  using  a  BATCH  submittal 
procedure,  RPLTQWM.  The  following  lines  must 
be  edited  for  the  specification  of  input  and  output 
file  names. 

GET,  PLTWQM,  input  file  name 

LGO,  input  file  name,  OPL 

1 ,  CALL,  SAVER(I=PLFILE,  J=plot  output  file 

name) 
REPLACE,  OPL=message  file  name 

The  input  file  is  a  monthly  data  table  (TAPE7 
output)  from  QWMON. 

The  procedure  is  run  by  the  interactive  commands: 

GET,  RPLTQWM 

(edit) 

SUB 

The  plot  output  contains  data  used  to  generate 
plots  using  a  DISSPLA  post-processor  procedure, 
such  as  COMP80  or  TEKTRON.  The  message  file 
contains  execution  information  generated  by  the 
DISSPLA  graphics  package,  run  cost  data,  and  a 
complete  dayfile.  It  should  be  checked  prior  to 
generating  hardcopy  plots. 
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Mission  of  the  Bureau  of  Reclamation 

The  Bureau  of  Reclamation  of  the  U.S.  Department  of  the  Interior  is 
responsible  for  the  development  and  conservation  of  the  Nation's 
water  resources  in  the  Western  United  States. 

The  Bureau's  original  purpose  "to  proviae  for  the  reclamation  of  arid 
and  semiarid  lands  in  the  West"  today  covers  a  wide  range  of  interre- 
lated functions.  These  include  providing  municipal  and  industrial  water 
supplies;  hydroelectric  power  generation;  irrigation  water  for  agricul- 
ture; water  quality  improvement;  flood  control;  river  navigation;  river 
regulation  and  control;  fish  and  wildlife  enhancement;  outdoor  recrea- 
tion; and  research  on  water-related  design,  construction,  materials, 
atmospheric  management,  and  wind  and  solar  power. 

Bureau  programs  most  frequently  are  the  result  of  close  cooperation 
with  the  U.S.  Congress,  other  Federal  agencies.  States,  local  govern- 
ments, academic  institutions,  water-user  organizations,  and  other 
concerned  groups. 


A  free  pamphlet  is  available  from  the  Bureau  entitled  "Publications 
for  Sale."  It  describes  some  of  the  technical  publications  currently 
available,  their  cost,  and  how  to  order  them.  The  pamphlet  can  be 
obtained  upon  request  from  the  Bureau  of  Reclamation,  Attn  D-922, 
P  0  Box  25007,  Denver  Federal  Center,  Denver  CO  80225-0007. 
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